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ABSTRACT 


Procedures to Predict Vertical Differential Soil Movement 
For Expansive Soils. (December 1997) 
Donald David Naiser, Jr., B.S., Texas A&M University 
Chair of Advisory Committee: Dr. Robert L. Lytton 


Damage to lightly loaded structures, paving and service piping in areas of 
expansive clay soils has occurred throughout the world. The cause of this damage has been 
the inability to accurately model expansive soil movement so that foundations are 
adequately designed to withstand the movement. The amount and rate of differential soil 
movement for expansive soils is due to a combination of soil characteristics, namely: 
suction compression index, unsaturated permeability, and diffusivity. Currently, 
geotechnical engineers run tests to measure the soil properties required to estimate 
differential soil movements. However, there seems to be apprehension toward attempting 
these soil movement calculations due to the perceived complexity of the calculations or a 
simple lack of understanding of the theory. The procedures delineating the step by step 
process used to calculate suction profiles and volume strains of expansive soils 1s 
presented. These procedures include the methodology to predict soil heave and shrink 
underneath shallow foundations which generate maximum center lift and maximum edge 
lift slab distortion modes. The main contributions of this research are: equations and 
procedures to calculate the equilibrium suction profile and depth to constant suction for a 
particular soil profile and location, equations to calculate the horizontal velocity flow of 
water in unsaturated soils, the methodology to predict differential soil movement shortly 
after a slab has been constructed and before the soil under the slab has reached an 
equilibrium moisture content, and the procedures to apply differential soil movement 
theory to soil profiles with multiple layers and moisture effect cases to be used for shallow 


foundation design. 
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Mitchell’s unsaturated permeability for the i soil increment (pg 11) 
Depth to the bottom of layers (pg 48) 

Liquid limit for layer s, in percent (pg 48) 

Liquid limit for the i" soil increment, in percent (pg 21) 

Number of layers in the b" boring sample (pg 48) 

Number of suction measurements in a sample (pg 48) 

Ratio of suction values between the top and bottom layers within the i" 
increment (pg 56) 

Plasticity Index for the i" soil increment, in percent (pg 20) 

Plastic limit for layer s, in percent (pg 48) 

Plastic limit for the i" soil increment, in percent (pg 20) 


Mitchell’s unsaturated permeability for the i" soil increment (pg 25) 
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Universal gas constant (pg 13) 

Slope of the line between suction values h, and h,,, (Appendix A) 
Slope of the suction vs gravimetric water content wet line for the i” 
increment (pg 14) 

Conversion of TMI, T = TMI + 60 (pg 30) 

Local annual Thornthwaite Moisture Index, (pg 30) 

An assumed value of TMI at the moisture depth d, (pg 30) 

Absolute temperature, degrees kelvin (pg. 13) 

Suction amplitude at depth L, (pg 57) 

Suction amplitude at the top of the incremental layer being considered 
(pg 36) 

Horizontal velocity of soil water at depth z, (pg 40) 

X coordinate for the b" boring sample (pg 48) 

Y coordinate for the b" boring sample (pg 48) 

Term used to damp the under slab suction profile toward the equilibrium 
suction profile (Appendix A-111) 

The sequential number describing a boring sample (pg 48) 

Suction measurement counter (pg 48) 

Mean annual moisture depth, d,, (pg 30) 

Available moisture depth of the soil mass (pg 30) 

Available moisture depth of the soil mass associated with the dry limit 
(pg 33) 

Size of horizontal barrier (pg 48) 

Mean annual moisture depth (pg 33) 

Distance between a moisture effect source and edge of barrier/slab (pg 46) 
Depth of vertical barrier (pg 48) 

Available moisture depth of the soil mass associated with the wet limit 
(pg 33) 

Depth of the root or flower bed zone (pg 46) 
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Moisture depth used to arrive at T, (pg 30) 

Controlling center lift, edge moisture variation distance (pg 52) 

Center lift, edge moisture variation distance for the i" soil increment (pg 27) 
Controlling edge lift, edge moisture variation distance (pg 52) 

Edge lift, edge moisture variation distance for the i" soil increment (pg 27) 


Edge moisture variation distance for return period r (pg 27) 
Edge moisture variation distance for a 10 year return period (pg 27) 
Edge moisture variation distance for a 50 year return period (pg 27) 


Soil void ratio (pg 25) 

A functional relationship between moisture depth d and T (pg 30) 
Cracked fabric factor for the i" soil increment (pg 20) 

Conversion from gram mass to gram force (pg 13) 

total suction, tension (usual units, pF) (pg. 13) 

c suction measurement in a sample (usual units, pF) (pg 48) 

Soil suction value at the time of construction (pF units) (pg 44) 

Soil suction value at the time of construction for the i" soil increment (pF 
units) (pg 45) 

Dry suction limit for soils in the field, bare soil at the surface usually 
has a value of 6.0pF and soil with vegetation at the surface usually 

has a value of 4.5pF (pF units) (pg 14) 

Soil suction profile at the edge of the slab at depth z, (pF units) (pg 63) 
Vertical soil suction profile at the entry point between two locations 
separated by a horizontal distance, when calculating the horizontal water 
flow (pF units) (pg 41) 

Vertical soil suction profile at the exit point between two locations 
separated by a horizontal distance, when calculating the horizontal water 
flow (pF units) (pg 41) 

Final suction value at depth z, (pF units) (pg 19) 
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Suction profile with respect to time and depth for a flower bed at the 
surface (pF units) (Appendix A-110) 

Field capacity suction value of soil (pg 14) 

Suction profile with respect to time and depth for grass at the surface 
(pF units) (Appendix A-110) 

Initial suction value at depth z, (pF units) (pg 19) 

Equilibrium suction value of the soil mass (pF units) (pg 30) 
Equilibrium suction value at depth z, (pF units) (pg 30) 

Dry limit suction for soils in the field, bare soil at the surface usually 
has a value of 6.0pF and soil with vegetation at the surface usually 


has a value of 4.5pF (pg 30) 


Dnax(Z)/hy,(Z,) Dry limit suction value at depth z, (pF units) (pg 30) 
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Wet limit suction of soils in the field, usually 2.5pF (pF units) (pg 30) 
Wet limit suction value at depth z, (pF units) (pg 30) 

Suction profile with respect to time and depth for soil underneath a slab 
(pF units) (Appendix A-110) 

Suction profile with respect to depth of under slab soil at time of 
construction (pF units) (pg 42) 

Suction profile with respect to depth of under slab soil at time the soil 
volume change is sought (pF units) (Appendix A-111) 

Suction profile with respect to time and depth for bare soil at the surface 
(pF units) (Appendix A-110) 

Suction profile with respect to time and depth for a tree at the surface 
(pF units) (Appendix A-110) 

Wet limit suction for soils in the field, usually 2.5pF (pg 14) 

Suction profile with respect to time and depth (pF units) (pg 30) 

Soil suction midway between the suction value at time of construction and 
the equilibrium suction value (pF units) (pg 44) 

Soil increment, 5 cm. (Appendix A) 
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Incremental layer counter (pg 19) 

Lateral earth pressure coefficient (pg 19) 

Gram molecular weight of water, 18.02 gm/mole (pg 13) 

Local annual weather cycle (pg 36) 

Return period (in years) (pg 27) 

Dummy variable (pg 34) 

Soil layer counter (pg 48) 

Time in seconds (pg 38) 

The month the slab was constructed (pg 45) 

Time elapsed between the driest month and the month of construction 
(pg 66) 

Driest month of the year (pg 49) 

Dummy time variable (Appendix A-107) 

Dummy time variable (Appendix A-107) 

Dummy time variable (Appendix A-107) 

The month the soil volume change is sought (pg 49) 

Time elapsed between the month of construction and the month soil volume 
change is sought (pg 66) 

Time elapsed between the driest month of the year and the month soil 
volume change is sought (pg 66) 


Time after construction when the value h,(z,) is reached (pg 44) 
Partial pressure of pore water vapor (pg x1) 


Saturation pressure of water vapor over a flat surface of pure water at the 
same temperature (pg x1) 

Horizontal distance (pg 40) 

The characteristic soil depth over which the available moisture depth is 
removed to achieve the driest soil state (pg 32) 


Depth to the i" soil increment (pg 19) 
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Depth to constant suction (pg | 1) 

Scores computed for return period r for the Gumbel cumulative probability 
distribution curve (pg 27) 

Depth to the top of the incremental soil layer being considered 
(Appendix A) 

Scores computed for a 10 year return period for the Gumbel cumulative 
probability distribution curve (pg 27) 

Scores computed for a 50 year return period for the Gumbel cumulative 
probability distribution curve (pg 27) 

Mitchell’s diffusion coefficient for the i” soil increment (pg 11) 

Shape factor in Gumbel’s cumulative probability equation (pg 28) 


An exponent from Mitchell’s equation set equal to 1 (pg 44) 


Vertical strain for the i” soil increment (pg 20) 
Volumetric strain due to suction for the i" soil increment (Appendix A-103) 


Volumetric strain correction due to overburden for the i" soil increment 


(Appendix A-103) 


Volumetric strain due to suction and overburden for the i" soil increment 


(pg 19) 

Total heave or shrinkage in the soil mass (pg 20) 

Vertical heave or shrinkage in the soil mass at the edge of the slab (pg 65) 
Vertical heave or shrinkage in the soil mass at the edge moisture variation 
distance measured to the center of the slab from the edge of the slab (pg 65) 
The difference in vertical movement between the edge of the slab and a 
vertical movement at the edge moisture variation distance measured to the 
center of the slab from the edge of the slab (pg 65) 

Regression coefficient (pg 30) 
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Dry unit weight for layer s (Appendix A-84) 

Dry unit weight for the soil in the i" increment (pg 25) 

Matrix suction compression index for the i" soil increment (pg 11) 
Matrix suction compression for a soil consisting of 100% of a particular 
type of clay mineral (pg 52) 

Osmotic suction compression index for the i" soil increment (pg 19) 
Mean principal stress compression index for the i soil increment (pg 19) 
Total unit weight (pg 19) 

Unit weight of water (pg 25) 

Matrix suction (pg 13) 

Percent passing 2 micron for layer s (Appendix A-84) 

Percent passing 2 micron for the i" soil increment (pg 21) 

Percent passing number 200 sieve for layer s (Appendix A-84) 

Percent passing the number 200 sieve for the i" soil increment (pg 21) 
Osmotic suction (pg 13) 

Final osmotic suction value a depth z; (pg 19) 

Initial osmotic suction value a depth z, (pg 19) 

Volumetric water content at depth z,, corresponding to hy, (pg 32) 
Volumetric water content at depth z,, corresponding to h,, (pg 32) 
Residual volumetric water content for the soil at depth z,, (pg 33) 
Saturated volumetric water content for the soil at depth z,, (pg 33) 
Volumetric water content at depth z,, corresponding to h,, (pg 32) 
Shape factor in Gumbel’s probability equation (pg 28) 

Applied octahedral normal stress at depth z, (pg 19) 

Stress level below which there is no suppression of suction volume strain by 


overburden pressure (pg 19) 


Total suction (pg 13) 





Activity ratio 


Cation exchange 


Cation exchange activity 


Cation exchange capacity 


Center lift distortion mode 


Clay 


Depth of available moisture 


Depth to constant suction 
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LIST OF DEFINITIONS 


The ratio of the plasticity index to the percentage clay 
content (Skempton 1955). Percent clay is defined as the 
percentage less than 2m of the soil material passing the U.S. 
No. 200 sieve (Pearring 1968). 

The interchange between a cation in solution and 

another cation on a surface active material (Holt 1969). 

The ratio of the cation exchange capacity to the percentage 
clay content (Pearring 1968). 

The total amount of exchangeable cations that a soil 

is capable of absorbing, measured in milliequivalents per 100 
grams of soil (Mojekwu 1979). 

This is a slab-on-ground distortion mode resulting from the 
soil at the perimeter being more dry (perimeter soil may be 
shrinking) than the soil beneath the slab (soil beneath the 
slab is either at equilibrium, shrinking at a rate slower than 
the exterior soil, or heaving). Refer to Fig. 2 (Castleberry 
1974 and Mitchell 1980). 

The soil particles smaller than 2. which are derived from 
the chemical decomposition of rock. It consists, chiefly, of 
clay minerals but small amounts of quartz, feldspar, organic 
matter, soluble salts and amorphous materials are also 
present (Lambe and Whitman 1969). 

The maximum depth of moisture available for use by 
transpiring vegetation, which is stored within the soil zone 
up to the depth to constant suction (Ga.y 1994). 


The depth in a soil profile to which there is no longer a 





Edge lift distortion mode 


Edge moisture variation 


distance 


Equilibrium suction profile 


Horizontal barrier 


Liquid limit 


Mean principal stress 


X1X 


significant seasonal suction change (the seasonal change in 
suction is less than 0.2 pF (Lytton 1997)). 

This is a distortion mode resulting form the soil at 

the perimeter being more wet (perimeter soil may be 
heaving) than the soil beneath the slab (soil beneath the slab 
is either at equilibrium, heaving at a rate slower than the 
exterior soil, or shrinking). Refer to Fig. 3 (Mitchell 1980). 
For soil under an impervious membrane, the edge 

moisture variation distance is the distance from the edge of 
the barrier to the point where soil suction does not vary with 
season moisture variation (Jayatilaka et al. 1992). 

A suction profile generated by decreasing the equilibrium 
suction, h,, a centimeter in suction for every centimeter 
increase in depth (McQueen and Miller 1968). 

Any impervious barrier that is attached and sealed to 

the edge of the foundation and extends some distance, d,,,, 
horizontally from the edge of the slab. A sidewalk attached 
and sealed to the edge of a foundation is considered a 
horizontal barrier. Refer to the sidewalk shown in Appendix 
C-132. 

The water content, expressed as a percentage of the weight 
of the dry soil, at the boundary between the liquid and plastic 
states. The water content at this boundary is defined as the 
water content at which a groove, 1mm wide, cut in the soil 
sample with a standard grooving tool, closes for a length of 
’2 inch when a brass dish containing the soil is dropped 25 
times, at a rate of 2 drops per second, through a distance of 1 
cm on a standard hard rubber base (Mojekwu 1979). 


The slope of the volume-total stress curve in the plane of 





compression index 


Plastic limit 


Plasticity index 


Post construction 


Post equilibrium 


Soil suction 


Suction amplitude 


Suction compression index 


Suction profile 


XX 


zero suction (Lytton 1994). 

The water content, expressed as a percentage of the dry 
weight of the soil at which the soil becomes crumbly and 
ceases to be plastic. The water content is defined as the 
water content at which the soil just begins to break apart and 
crumble when rolled, by hand, into threads one-eighth of an 
inch in diameter (Mojekwu 1979). 

The difference between the greater moisture content 

of the liquid limit and the less moisture content of the plastic 
limit. It is defined as the range in water contents through 
which the soil remains in the plastic state (Mojekwu 1979). 
This refers to the transient case for calculating volume 
strain of expansive soils immediately after construction, 
before the soil under the slab has reached an equilibrium 
moisture content. 

This refers to the case for calculating volume strain of 
expansive soils after the soil under the slab has reached an 
equilibrium moisture content. 

the negative gage pressure relative to the external gas 
pressure on the soil water to which a pool of pure water 
must be subjected in order to be in equilibrium through a 
semi-permeable membrane with the soil water. 

The absolute value of the range between the maximum and 
the minimum suction values being considered in soil 
suction profile, (typically, U, =h,,,,(Z;) - Mmin(Z;,) OF 

U,= h(z) - h,,(2)). 

The slope of the volume-total suction curve in the plane of 
zero suction (Lytton 1994). 


When the state of suction is measured at intervals of depth 





ere 


down the soil profile, the resulting relation between suction 
and depth is termed the suction profile (Mitchell 1981). 

Thornthwaite Moisture Index TMI is defined as a climatic index that measures the water 
balance of an area to determine if there is a deficiency or 
excess to ascertain if the climate should be classified as arid 
or humid (Thornthwaite 1948). 

Vertical barrier Any impervious barrier that extends some distance, d,,,,. 
vertically downward from the foundation edge. 

Vertical differential The total amount of vertical movement between two 


Soil Movement distinct locations separated by a characteristic distance. 





INTRODUCTION 


When a slab barrier is placed on an expansive soil, the immediate effect is that 
evaporation and precipitation cannot effect the soil beneath the slab. During subsequent 
weather events occurring to the surrounding soil, moisture differentials develop between 
the uncovered soil and the soil beneath the slab barrier. These differentials are largest 
when the climate is described by a distinct wet period followed by a hot dry period. Unless 
a foundation system is installed to prevent or reduce moisture movement, moisture 
redistribution occurs within the soil mass to reduce these differences in moisture content 
(Holland and Lawrance 1980). Damage caused by this type of moisture differential is 
considerable and warrants investigation in developing a more scientific and accurate 
methodology to address these problems. 

The ability to predict differential soil movement is critical for engineers to 
formulate a picture of what the soil/structure interface will be. This information is 
necessary to adequately design foundations or other ground supported structures. In slab- 
on-ground design, the principal interest is in making an accurate estimate of the range of 
vertical movement that must be sustained by the foundation. The procedures contained 
herein are methods for defining envelopes of maximum heave and shrinkage so that 
adequate slab-on-ground designs are made. For highway and airport pavements, canals, 
and pipelines, the wave spectrum of differential movements versus wave lengths are the 
desirable design characteristics. Structural floors suspended above expansive clays must 
be provided with a gap that exceeds the total expected heave. Drilled piers (or shafts) must 
be designed to resist simultaneously a vertical movement profile and a horizontal pressure 
profile, both of which change with wetting and drying conditions. Retaining structures, 
basement walls, rip rap, and canal linings must be designed to withstand lateral 
movements. Finally, all foundations must be designed against the time-dependent vertical 


and horizontal curvature that is generated by down hill creep (Lytton 1994). 


This thesis follows the style of the Journal of Geotechnical Engineering. 





Each type of soil movement listed above and their effects on different types of 
ground structures is of sufficient importance and complexity to warrant research of its own. 
Addressed herein are the effects of vertical differential soil movement on commercial and 
residential types of slabs-on-ground. Vertical differential soil movement in these cases is 
defined as the total amount of vertical movement between two distinct locations separated 
by a characteristic distance. The characteristic distance is described as the edge moisture 
variation distance plus the distance a moisture effect source is away from the slab edge. 
The model used to predict differential movement of expansive soils is based on constitutive 
soil relationships derived primarily from works developed by C. W. Thornthwaite (1948), 
W.R. Gardner (1958), R. L. Lytton (1973, 1977, 1994, 1995, 1997), J. R. Pearring (1968), 
J. H. Holt (1969), E. C. Mojekwu (1979), J. L. Nieber (1981), R. G. McKeen (1981, 1992), 
P. W. Mitchell (1980, 1984), D. A. Gay (1994). 

It has been well documented that the magnitude at which soil shrinks and swells is 
not exclusively based on the soils characteristics. The magnitude of differential soil 
movement in expansive soils is dominated by the following variables: soil characteristics, 
pre-construction and post construction vegetation, slope of perimeter surface grade, shape 
and age of slab and depth of soil moisture active zone (Castleberry 1974, Mathewson et al. 
1975, Mathewson et al. 1980). 

The methodology used to predict vertical differential soil movement involves using 
soil suction profiles under the foundation and at the edge of the foundation. These suction 
profiles are developed using the following variables: depth, time, local surface annual 
weather and vegetation conditions, and soil characteristics such as the suction compression 
index, unsaturated permeability, and soil diffusivity. Soil suction profiles combined with 
overburden stresses are used in volume strain equations developed by Lytton (1994) to 
predict the volumetric strain for an incremental layer of soil. The incremental layer 
movements are summed up and total vertical differential soil movement is predicted for a 
particular column of soil. These procedures are detailed in the chapter on “The Procedures, 
in a Software Design Format”. These procedures address foundation perimeter moisture 


effect cases that are either very common with many foundation problems or have a 





significant effect on vertical differential soil movement. Specifically, the foundation 
perimeter moisture effect cases that are addressed include: (1) bare soil at the surface, (2) 
grass at the surface, (3) flower bed at the surface with a defined moisture effect zone depth, 
(4) trees at the surface with a defined root zone depth, and (5) moisture effect zones that 
have been measured in the form of suction profiles. Included are methods for predicting the 
effects that vertical or horizontal barriers have on vertical differential soil movement. 
Additionally, methods to predict differential soil movement shortly after construction and 
before the soil has reached an equilibrium moisture content under the foundation are the 
highlights of this research and are presented. 

In “The Problem” chapter, the theory and methodology used to predict differential 
soil movement is defined and described in order by which the theory was developed. 
These bits and pieces of theory are applied in step by step processes described in the 
procedures sections of this document to address all of the different moisture effect cases, 


with and without barriers imposed. 





THE PROBLEM 


The Magnitude and Description of the Problem 


Since early 1950, damage to lightly loaded structures, paving and service piping in 
areas of expansive clay soils has occurred throughout the world. Estimates show that one 
in five Americans are affected by this type of damage while only one in ten are affected by 
major floods. Structural damage caused by expansive clay soils in the United States 
annually exceeds that caused by earthquakes, hurricanes, and floods combined (Jones and 
Holtz 1973.) Damage to residential and commercial structures caused by soil movement is 
estimated at 660 million dollars, almost one-third of the total monetary damage received by 
all structures as a result of soil movement. On a state-by-state basis, California and Texas 
account for 35% of all expansive soil losses in the United States. This becomes obvious 
from a map of the location of much of the expansive soils in the United States, see Fig. 1. 


Much of this soil is located in highly populated regions of both California and Texas. 
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Fig. 1. United States Map of Expansive Soils After Wiggins (1976) 





When it is realized that most residences and light commercial buildings in these two states, 
which together had 430,000 new home starts in 1978, are constructed with slab-on-ground 
foundations, the need to develop a theoretical design approach based on physical and 
natural law for slab-on-ground design in expansive soil areas becomes evident. 

The problem with slabs-on-ground built on expansive soils becomes manifested as 
excessive structural distortions in response to the subsoil undergoing volume changes 
induced by a change in subsoil moisture. If foundations are not of a stiffness to resist or 
span the gaps caused by differential soil movements, and if the superstructure is not of a 
suitable flexibility, the foundation distortions resulting from the soil movements can be of 
a magnitude to cause excessive cracking of walls and distortion of internal fittings. This 
type of differential shrink and heave of soils give rise to two distinct distortion modes 
which cause damage to the superstructure. Fig 2 shows typical superstructure damage that 
would occur in the case of center lift. Center lift is a distortion mode resulting from the 
soil at the perimeter being more dry (perimeter soil may be shrinking) than the soil beneath 
the slab (soil beneath the slab is either at equilibrium, shrinking at a rate slower than the 
exterior soil, or heaving.) Fig 3 shows typical superstructure damage that would occur in 
the case of edge lift. Edge lift 1s a distortion mode resulting from the soil at the perimeter 
being more wet (perimeter soil may be heaving) than the soil beneath the slab (soil beneath 
the slab is either at equilibrium, heaving at a rate slower than the exterior soil, or 
shrinking.) The engineering profession lacks a universal method of analysis within a 
theoretical framework that utilizes fundamental properties of soil and footing behavior to 
enable a footing design that is satisfactorily accurate without resorting to empirical 
methods (Mitchell 1980). Since the late 1970’s strides have been made by the Post 
Tensioning Institute (PTI) to adopt a standard method for designing ground structures in 
expansive soils. The Post Tensioned Institute has adopted the methodology discussed 
herein and stipulates that post tensioned slab designs shall provide necessary rigidity to 
avoid damage to the buildings or houses when the supporting soil expands or contracts 
(PTI 1996). However, a theoretical approach to predict vertical differential soil movement, 


as presented herein, has not been widely adopted by today’s practicing engineers. 
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Fig. 2. (Mitchell 1980) Center Lift Distortion Mode Fig. 3. (Mitchell 1980) Edge Lift Distortion Mode 


Expansive Soil, on a Molecular Level 


The Building Blocks of Expansive Soil 


Why are “expansive soils” so expansive? The answer can be found in the material 
on a molecular level. The magnitude of the surface area per mass and the molecular 
composition of expansive soils plays a major role in the soil’s ability to shrink and swell to 
such an impressive degree. “Expansive soils” are typically the fine grained soils as 
classified in the unified classification system. Fifty percent or more of a soil, by weight, in 
a sample that passes the number 200 sieve is classified as a fine grained soil. Fine grained 
soils are further broken down into categories of silts, clays and organics having either a 
high or low plasticity. The ability to break down and classify a soil according to the unified 
soil classification system provides an engineer with a wealth of information about the 
behavioral characteristics and how a particular soil may effect design. An in depth 
explanation of soil characteristics and the unified soil classification system can be found in 
Lambe and Whitman (1969). 

The following clay structural unit explanation is idealized and simplified and are 
presented to provide a picture of the basic behavioral characteristics that give rise to clay’s 
ability to undergo large changes in volume. The structural units that make up a clay 


particle are silicate sheets. Primarily, clays are made of sandwiches of: silica sheets, 





gibbsite sheets, and brucite sheets. These sheets are combined in two distinct groups: 
two-layer sheet minerals, and three layer sheet minerals. The most important and most 
common two-layer sheet mineral 1s the kaolinite clay particle, see Fig. 4. Kaolinite is a 
two-layer mineral which consists of a silicate sheet bonded to a gibbsite sheet by 
hydrogen bonding and secondary valence forces. Another example of a two layer sheet 


mineral is serpentine which is formed by brucite bonding to silica, see Fig. 5. 
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Fig. 4. (Lambe and Whitman 1969) The Structure of Kaolinite 
(a) Atomic Structure (b) Symbolic Structure 


In the actual formation of these silicate sheets there may be some isomorphous 


substitution. Isomorphous substitution is when one kind of atom 1s substituted for 
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Fig. 5. (Lambe and Whitman 1969) The Structure of Serpentine 
(a) Atomic Structure (b) Symbolic Structure 


another atom that is very similar in structure and size. In other words aluminum may be 
substituted for silicon if there is a shortage of silicon and an abundance of aluminum 
during the formation stage of the clay mineral. This sort of substitution may give the 
mineral a net charge deficiency and may cause slight distortion in the crystal lattice due 
to the substituted atom’s size difference. Isomorphous substitution that occurs giving 
the clay mineral a large negative charge will give the mineral the capacity to attract 
polar charged water molecules. A double layer of molecules, such as water molecules, 
will try to form around a clay mineral to neutralize the negatively charged clay particle. 
This great affinity for water due to isomorphous substitutions provides an explanation 


as to why clay minerals, like montmorillonite, can have large volume changes at high 





volumetric water contents. In Fig. 6 the clay particles are shown with the fully 


developed double layers they would have in pure water. 
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Fig. 6. (Lambe and Whitman 1969) Soil Particles With Water and Ions 
(a) Sodium Montmorillonite (b) Sodium Kaolinite 

In kaolinite there is a very small amount of isomorphous substitution. However 
in the three layer sheet clay minerals there is much more isomorphous substitution. 
Three layer sheets are formed by placing a silica on the top of, and on the bottom of, 
either a gibbsite or brucite sheet. Fig. 7 shows the mineral pyrophyllite made of a 
gibbsite sheet sandwiched between two silica sheets. Fig. 8 shows the structure of the 
mineral muscovite, which is similar to pyrophyllite except that there has been 
isomorphous substitution of aluminum for silicon in muscovite. The net charge created 
by this substitution is balanced by potassium ions which serve to link the three layer 
sandwiches together. The most common three-layer minerals in soil are 
montmorillonite and illite. Montmorillonite is similar to pyrophyllite with the 
exception that there has been isomorphous substitution of magnesium for aluminum in 
the gibbsite sheets. Due to the isomorphous substitution the illite mineral has a 
moderately high capacity for volume change and montmorillonite has a very high 


capacity for volume change (Lambe and Whitman 1969). 
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Fig. 7. (Lambe and Whitman 1969) The Structure of Pyrophyllite 


(a) Atomic Structure (b) Symbolic Structure 
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Fig. 8. (Lambe and Whitman 1969) The Structure of Muscovite 


(a) Atomic Structure (b) Symbolic Structure 





Gilgai 

Clay soil left to interact with surface climate conditions will develop an irregular 
surface microrelief pattern called gilgai, which in colloquial terms is referred to as “hog 
wallows.” This phenomenon is recognizable in areas where highly expansive soils are 
known to exist undisturbed for extended periods of time. These gilgai may be 
experienced while taking a very undulating ride: up the fairway of an old golf course, 
during farming evolutions in a native hay field, or down an old country road. As 
illustrated in Fig. 9 from Spotts (1974), the pattern develops in seven stages. First, the 
sedimentary clay is deposited. Following deposition, drying occurs, producing surface 
shrinkage cracks. Further desiccation results in the formation of major cracks 
penetrating far into the clay. The extent of penetration and spacing are determined by a 
number of factors of which the degree and rate of drying and the clay type are the most 
significant. Subsequent development of microrelief occurs as the clay is wetted and 
dried by cyclic variation of the climate. Because water penetrates deeper at the major 
cracks, greater swell occurs. As the next cycle of drying proceeds, the higher areas dry 
first, causing the major cracks to reappear in approximately the same locations. Thus, 
over a period of several years, an undulating surface pattern called gilgai develops on 
the soil surface (McKeen 1981). This phenomenon is important to realize since it is a 
natural occurrence of differential soil movement that must be understood by engineers 
so that appropriate models of this type of soil behavior can be properly applied to 
design. There are important parameters to consider when addressing volume strain that 
causes the gilgai phenomenon, such as: (1) The depth at which soil interacts with the 
climate, known as the depth to constant suction, z,,. (2) The climate, or the annual 
wetting and drying cycles at the surface. (3) The type of surface conditions, such as; 
vegetation growth, bare soil, slope, or barriers of some type. (4) The soil characteristics, 


namely: suction compression index, y,(z;), unsaturated permeability, K,(z,)|h,|, and soil 





hie 





diffusivity, a(z,). The roles that each of these parameters play in predicting differential 


soil movement are presented. 
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Fig. 9. (Spotts 1974) Gilgai Stages 





Suction, Matrix and Osmotic 


Soil suction is the measure of a soil’s affinity for water. Soil suction is an 
energy: water always moves down an energy gradient and it takes energy to cause 
volume strain in any soil. Soil scientists and arboriculture scientists refer to soil suction 
in terms of water potential. As defined by the 1960 Moisture Equilibria Symposium, 
soil suction is the negative gage pressure relative to the external gas pressure on the soil 
water to which a pool of pure water must be subjected in order to be in equilibrium 
through a semi-permeable membrane with the soil water. Suction is the term used 
principally by engineers for the thermodynamic quantity, Gibb’s free energy which is 
inherently negative, as seen in (1) below, and generates tension in the pore water 
stretching between soil particles (Lytton 1994). Soil suction is equal to 0.0 when the 

= aa inf) (1) 
m-g 100 
relative humidity is equal to 100%. A relative humidity value less than 100% in a soil 


would indicate the presence of suction in the soil. 


Components of Soil Suction 


Soil suction is quantified in terms of the relative humidity (1) and is commonly 
called “total suction.” “Total suction,” y , has two components, namely: matrix suction 
(1), which is due to the attraction of water to the soil particle surfaces, and osmotic 
suction (7), which is due to dissolved salts or other solutes in the pore water (2). 

YW=UrtT. (2) 
Matrix suction will have the most effect on the volumetric changes in soil due to the fact 
that this component varies with the changes in the moisture environment. Osmotic 
suction will not change unless there is a change in the salt concentration within a given 


soil layer. 


Typical Suction Levels in the Field 


Fig. 10 illustrates the suction-vs-volumetric water content curve for a natural soil under 





wetting and drying conditions. A common measure of suction is the pF-scale, in which 
pF is defined as (3), where p refers to the logarithmic value analagous to pH, and the F 
refers to free energy. 

Suction in pF = log,| h| (3) 
Where | h| equals the magnitude of suction in centimeters of water, a positive value. 


The various suction levels corresponding to field cases are listed below: 
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Fig. 10. (Lytton 1994) Suction vs Volumetric Water Content Curve 


Field capacity (h,, = 2.0 pF) 

Wet Limit for clays (h,,= 2.5 pF) 

Plastic limit (3.5 pF) 

Wilting point of plants (hy, = 4.5 pF for vegetation) 

Tensile strength of confined water (5.3 pF) 

Air dry at 50% relative humidity (h,,, = 6.0 pF for bare soil) 
Oven dry (7.0 pF) 





Field capacity is an agricultural engineering term used to define the maximum amount 
of water that a soil can hold in the field. Fig. 10 shows field capacity to be a value of 
2.0 pF, but actual field values never quite drop below a value of 2.5 pF for clays (Lytton 
1995 and Mitchell 1980). Soil suction is a term that measures energy potential and is 
directly related to the soil’s potential for volume change. This information is referred to 


repeatedly throughout this document. 





THE THEORY OF DIFFERENTIAL SOIL MOVEMENT 


Some notable breakthroughs have been made pertaining to the topic of 
predicting differential soil movement and are described in this chapter. The real 
challenge in geotechnical engineering, especially in the area of expansive soils, is to 
develop a consistent body of theoretical mechanics which adequately explains the 
behavior of the soil mass and to develop methods of measuring the relevant material 
properties critical for making appropriate design decisions. 

In the 1970’s the first attempt was made to design slabs-on-ground using a 
scientifically based theoretical approach. In 1973, Lytton and Woodburn prepared a 
report which chronicles a successful school building foundation design. The school 
foundation design procedure considered differentially expanding and shrinking clay due 
to the changing suction of the soil beneath the foundation. The suction profiles of the 
soil were predicted using a rational method which relied upon suction profiles measured 
on samples taken on the site at different times in the year. These suction values were 
used to predict the differential soil movement, and the foundations were designed based 
on the distortion modes caused by the magnitude and location of the expansive soil 
movements. Lytton noted that this differentially expanding and shrinking effect caused 
“mounds” to form beneath the slab. The foundation interacts with these mounds 
pressing down on the high spots and bridging the low spots. The foundation design 
procedures considers the limit cases in which the mounds produce two types of 
distortion patterns: edge lift pattern and center lift pattern, as shown in Fig. 11. These 
two distortion modes were used to design the slabs so that the loads of the 
superstructure would be carried across a span generated by edge lift, or carried in a 
cantilevered fashion out to the edge of the foundation in the case of center lift. 
Woodburn recorded the performance of the foundation systems after the designs were 
built and the school was put into service. Lytton and Woodburn reported post 
construction edge shrink values to be 2.75 inches which compared favorably with the 


predicted value of 3.1 inches. Just as Lytton predicted, the worst case differential 





movement occurred in the first dry season subsequent to the construction of the slab 
during the wettest time of the year. In this post construction case the soil was “capped” 


and the moisture content was preserved in a 





Edge Heave Center Heave 


Fig. 11. (Lytton 1973) Edge Heave and Center Heave Distortion Modes 


state of heave under the slab while the perimeter soil shrank as it dried. This post 
construction case produced the worst case center lift mode. Subsequent weather cycles 
did not produce such extreme cases of differential soil movements as compared to the 
post construction case. However, subsequent weather cycles did produce differential 
soil movements as predicted. The model to predict this type of post construction 
differential soil movement is presented in the Post Construction section of this chapter. 
This design and research venture documented by Lytton and Woodburn’s 1973 
paper marked the beginning of a slab-on-ground design method that has been adopted 
by the Post Tensioning Institute as the standard for design of post-tensioned slabs-on- 
ground on expansive soils (PTI 1980, 1996). However, the objective of this research is 
to consider the methodology behind predicting the differential soil movement used to 


design such slabs. 


Soil Volume Strain 


Movements in expansive soils are generated by changes of suction which are 


brought about by the entry or loss of moisture or change in concentration of dissolved 





salts in the pore water. The volume change that accompanies the change of suction (and 
water content) depends upon the total stress states that surround the soil. Within a soil 
mass, a decrease in the magnitude of suction results in an increase of water content. 
The volume of the soil also increases unless the surrounding pressure is sufficient to 
restrain the swelling. A conceptual graph of suction-vs-volume can be drawn using the 
relations of each to water content. This is illustrated in Fig. 12 on the plane 
corresponding to zero pressure. A similar graph can be drawn relating pressure (total 
stress) - versus - volume on the plane corresponding to zero suction. The simultaneous 
change of the magnitude of suction (decrease) and pressure (increase) results in a small 
change of volume, following the path from Point A to Point C on the pressure - suction 
- volume surface. The magnitude of suction decreases from Point A to Point B while 
the pressure increases from Point B to Point C. The volume change process can be 
viewed as the net result of two processes; 

a. Increase of volume from A to B at constant mechanical pressure or total 

stress. 


b. Decrease of volume from B to C at constant suction. 
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Fig. 12. (Lytton 1994) Mean Principal Stress-Volume-Matrix Suction Diagram 





For small increments of volume change on this surface, the volume strain, is linearly 
related to the logarithms of both pressure and |suction|. The soil volume strain equation 


which shows the relation between matrix suction, mean principal stress, and osmotic 


suction 1s (4), (Lytton 1977, 1994). 





AV Nez: az -(Z, 
ia) i. -14)logo| 2) 16 logy| 2222) 1,2 )logy| 222) (4) 


In the procedures described herein the assumption will be made that the osmotic 
component has no change. If the case exists that the osmotic component has no real 
change, which is the case most of the time, then that term within the volume strain 
equation can be neglected. The suction compression index , y,(z,), and the mean 
principal stress compression index, y,(z;), will be discussed in more detail in the section 
titled “The Suction/Mean Principal Stress Compression Indexes”. The initial and final 
suction profiles, h,(z,) and h,(z,), will be discussed in more detail in the section titled 
“Suction Profiles”. The initial mean principal stress, o,, 1s the stress level below which 
there is no suppression of volume strain by overburden pressure. This is the stress at 
which the curve in the pressure-volume plane in Fig 12 begins to depart from horizontal 
as pressure increases from zero and has been observed to be approximately 40 
centimeters of soil (Lytton 1994). The final stress at depth z; 6,(z;), 1s the mean 


principal stress at a depth z, at the time the volume strain is desired and is given by (5). 





1+2-ko 
6,(Z;) = 3 “Y1° 4; (3) 
where k, is: k, = 0.00 when the soil is badly cracked. 


k, = 0.33 when the soil is drying. 
k, = 0.67 when the soil 1s wetting. 


k, = 1.00 when the cracks are closed and the soil is swelling. 
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Once the volume strain for an incremental layer of soil is known the vertical strain for 
that layer is estimated using a cracked fabric factor, f, , as shown in (6). The cracked 


fabric factor 1s the percent soil volume strain that is directed vertically in a soil column. 


(FF) c 


McKeen (1981) back calculated vales of f, and found that f; = 0.5 when soils are drying 
and f,= 0.8 when soils are wetting. The total vertical differential movement in a soil 
mass, Az, 1s then found by summing the product of vertical strain and the incremental 


depth to which they apply as shown in (7). 


(AH) , 
Az = » a: -Increment (7) 


The volume strain theory previously described was developed by Lytton (1973, 1977, 
1994, and 1995). All symbols are defined in the List of Symbols. 


A Simplified Method for Identifying Predominant Clays 


There are three areas of soil characteristics that must be known before the 
volume strain theory can be applied, namely: the soil suction compression index, the 
soil suction profiles, and the mean principal stresses within the soil. The mean principal 
stress was described earlier and suction profiles will be discussed in a later section. 
However, the soil suction compression index, as described in the next section, depends 
on proper identification of the predominant clays in soils. A simplified method for 
identifying these predominant clays is discussed in this section. 

Before Lytton developed the methodology to predict vertical differential soil 
movement, Pearring (1968) and Holt (1969) had completed a correlation chart to aid in 
the identification of predominant clay mineral in a given soil. According to Pearring 
and Holt’s research, the two parameters used to classify the clay minerals are the cation 
exchange activity, CEAc(z,), and activity ratio, Ac(z,). Obtaining the measure of these 
two parameters requires the plasticity index, PI(z,), the cation exchange capacity, 


CEC(z,), and the percent of clay in the soil passing the number 200 sieve, Clay(z,), (8), 





2] 


(9), (10), (11). 











PI(z,) = LL(z,) - PL(q) (8) 
Clay(z,) = a (9) 
saan aa (0 
CEAc(z,) = —_ (11) 


All symbols are defined in the List of Symbols. Once the cation exchange activity and 
the activity ratio are known for a given clay mineral the type of clay is determined from 
Pie )3. 

The difficulty comes in arriving at the measure of cation exchange capacity. 
Until Pearring and Holt’s method of identifying predominant clays came to light the 
geotechnical engineer had to rely on experience or the more reliable, expensive, and 
arcane techniques used by clay mineralogist and soil physicists. However, Pearring 
and Holt’s method still required engineers to find the cation exchange capacity (CEC) 
of the soil, which is a soil property not typically evaluated in normal soils test 
laboratories. Mojekwu (1979) realized this problem and developed a simplified way to 
arrive at the CEC. His work involved data collected for a wide range and number of 
soil samples taken in the state of Texas. Among other data, Mojekwu collected and 
analyzed the Atterberg limits of these soils and found that very good correlations could 
be made that linked the plastic limit and the liquid limit to the CEC of a specific soil. 
The method he developed requires the simple Atterberg limits tests to find the plastic 
limit, PL, and the liquid limit, LL. Mojekwu ran a simple regression analysis, the 
dependent variable being CEC, the independent variables being PL and LL respectively, 
and produced very high R squared values of .9941 and .9942 for (12) and (13) 


respectively. 


CEC(z,) = PL(z,)"” (12) 
CEC@) =LL(z)" (13) 
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The larger of these two approximated values of CEC 1s used to find the CEAc, which is 
then used along with Ac to determine the predominant clay minerals in Fig. 13. These 
findings by Pearring, Holt, and Mojekwu paved the way for research by McKeen in 
developing a simplified method to identify the suction compression index, which is 


critical to predicting differential soil movements. 
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Fig. 13. (Pearring 1968 and Holt 1969) Chart to Identify Predominant Clays 


The Suction/Mean Principal Stress Compression Indexes 


Suction forces on a soil element are body forces rather than externally applied 
forces at the surface of the element. Because the suction stresses are dominating the 
effective stresses in expansive soils, volume change studies are made using total suction 
values. It is the slope of the volume change versus suction curve that quantifies soil 
response to moisture changes. Thus, the suction compression index, y,(z;), 1s defined as 


the slope of the volume-total suction curve, see Fig 14 (Lytton 1977). McKeen (1981) 





Ze 


analyzed soils data in which he measured volume strain corresponding to soil suction 
values to arrive at a corresponding suction compression index for particular types of 
expansive soils. In McKeen’s study, he measured the soil volume change using the 
conventional oedometer for cases when the soil was un-cracked and using the 
coefficient of linear extensibility (COLE) method when the clay was cracked. The 
measurements of soil suction were taken utilizing a variety of methods as well. The 
primary methods were the thermocouple psychrometer and filter paper methods. Using 


the basic concept 
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Fig. 14. (Lytton 1977) Pressure-Volume-Suction Surface 


described above, McKeen developed the chart method to determine y, by analyzing a 
large number of Soil Survey samples, adopting the methodology for determining 
predominant clays by Pearring (1968), Holt (1969) and Mojekwu (1979), and produced 


Fig 15 which gives values of Y,;9) without requiring suction tests. The value of Yy199 
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typically ranges from .033 to .220. This method requires only the activity ratio and the 
cation exchange activity of the soil to find the corresponding y, for particular soils. The 
chart method 1s used in the procedures to calculate y, in the software design format 
procedures. 

Another factor to consider in calculating soil volume strain is the effect of 
applied loads in suppressing the suction volume strain by overburden pressure. The 
mean principal stress compression index is required to complete the term to calculate 
the amount of suppression to the suction volume strain by overburden pressure. The 
mean principal stress compression index is calculated by (14) (Lytton 1977). 
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Fig. 15. (McKeen 1981) Chart to Determine the Suction Compression Index 
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The mean principal compression index is also calculated by the commonly used 
compression index, C,, by (15) (Lytton 1994). 
C 


wan l 
Yo Pes (15) 





Typical values of y, were calculated from results of swell tests and swell pressure tests 
reported by Vigayvergiya and Ghazzaly (1973) according to (14). The range of these 
calculated values of y, fell between 0.02 and 0.09 indicating that the ratio of y,/y, is 
typically at or slightly below 1.0. Therefore it is slightly conservative but convenient to 
assume y,, is equal to y, Until there is research that develops a more practical and cost 
effective method for obtaining more accurate values of y,, that coefficient will assume 
the value of y, for calculations of volume strain in the software design format 


procedures described in this document (Vigayvergiya and Ghazzaly 1973). 


Estimates of Unsaturated Soil Properties 


Before a clear explanation of the remaining sections is possible, estimates of 
unsaturated soil properties must be defined. Specifically, definitions are required for the 
following soil properties, namely: the slope of the suction versus volumetric water 
content wetting line shown in Fig. 10, S(z,), Mitchell’s diffusion coefficient, a(z,), and 
or P(z,). The slope of the suction versus 


Mitchell’s unsaturated permeability, kK, (z,)lh, 








volumetric water content wetting line is given by (16). Mitchell’s diffusion coefficient 
is given by (17). Finally, Mitchell’s unsaturated permeability is given by (18). The 
process flow sheet that describes procedures for calculating all soil properties required 


for the software design can be found in Appendix A-88. 








S(z,) = -20.29 + 0.1555-LL(z,)—0.117-PI(z,) + 0.0684- v(z, ) (16) 

a(Z, ) = 0.0029 — 0.000162 -S(z,) —0.0122-y,(z,) (17) 

K,(z,): h, Pa Aa came (18) 
S(z,) le 


The above equations were developed by Jayatilaka et al. (1992) and Lytton (1994). All 


symbols are defined in the list of symbols. 





26 


Edge Moisture Variation for Center Lift and Edge Lift Cases 


As discussed earlier, the effects of vertical differential soil movement on 
commercial and residential types of slabs-on-ground is dictated by the perimeter 
moisture effect cases and the total amount of vertical movement between two distinct 
locations separated by a characteristic distance. The characteristic distance is described 
as the edge moisture variation distance. At some time after construction, the soil under 
the slab will reach an equilibrium moisture content. However, the soil at the edge will 
vary in moisture content dependent on the local annual weather cycle. Additionally, the 
soil underneath the slab will vary in moisture content up to a distinct distance, the edge 
moisture variation distance, that is controlled by the permeability and diffusivity of the 
soil. This differential in moisture content of soil underneath the slab and the edge of the 
slab is what causes the center lift and edge lift distortion modes. As part of research 
conducted by Jayatilaka et al. (1992) a calibrated finite element program with coupled 
transient moisture flow and elasticity that had been used in the study of vertical 
moisture barriers provided an ideal means to study the edge moisture variation distance. 
The unsaturated soil properties discussed in the previous section in addition to the 
Thornthwaite Moisture Index were used to determine the relation of the edge moisture 
variation distance. 

The Thornthwaite Moisture Index (TMI), which measures the location’s 
potential for evapotranspiration, is defined as a climatic index that measures the water 
balance of an area to determine if there is a water deficiency or water excess. 

Both edge lift and center lift conditions were explored using several hundred 
runs with the program (Jayatilaka et al. 1992). Center lift conditions were simulated by 
a one year dry spell following a wet suction profile condition. Edge lift conditions 
were simulated by a one year wet spell following a dry suction profile condition. The 
edge moisture variation distance was considered to be that distance between the edge of 
the foundation and the point beneath the covered area where the suction changed no 


more that .2 pF during the entire period of simulation (Lytton 1994). 





2/ 


The dry and wet conditions used annual suction variation patterns that were 
appropriate for each of nine different climatic zones ranging from a Thornthwaite 
Moisture Index of -46.5 to +26.8, spanning the range found in Texas. The resulting 
edge moisture variation distances are shown in Figs. 16 and 17. No distance less that 
2.0 feet should be considered for design purposes (Lytton 1994). In the procedures for 
the software design, the edge moisture variation distances, e,,,(z,) and e,,.(Z,), are 
calculated using Figs. 16 and 17. 

In 1997 Lytton developed a method to determine the edge moisture variation 


distance for a particular design return period e,, . In designing slabs-on-ground 


recognition must be taken of the length of the time these structures must be in service, 
and of the severity of the weather patterns that may occur during the expected life of the 
structure. The return period of hydrologic events is appropriate to use in estimating the 
design criteria for slabs-on-ground on unsaturated soils. Determining the edge moisture 
variation distance for a selected return period involves applying an equation that relates 
the edge moisture variation values for the 10-year and 50-year return periods to the 


design edge moisture variation value for the selected return period, (19). 


eee. See (19) 
m; My Ms My 7. = 7 
50 10 


The edge moisture variation distance for a 10 year return period can be obtained from 
Fig. 18. These edge moisture variation distances were derived by back calculation from 
slabs which were performing successfully in San Antonio, Dallas, and Houston. None 
of the slabs were more than 10 years old at the time. It can be argued that the design 
values of edge moisture variation distance represent a 10 year return period, e,, . The 
edge moisture variation distance charts shown in Fig. 16 and Fig. 17 represent a 50 year 
return period, e,,,. The use of the Gumbel probability density function, which is 
commonly used to represent the probability of weather events, may be used to establish 


the risk level that is desired for design in accordance with the expected service life. 


The z, scores used in (19) are computed from the Gumbel cumulative probability 
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distribution curve, (20). The return period is represented by the variable r. p and f are 


z,=—_>—; (20) 


J) 


both shape factors and can be assumed as one. Solving for an edge moisture variation 
distance for a particular return period between 10 and 50 years is a simple two step 
process: (1) simply substitute the return periods of 10, 50 and the selected return 
period r into (20) and solve for Z,9, Z;) and z, respectively, and (2) Substitute values of 


e,,, from Fig. 18, e,,, from Fig. 16 or Fig. 17, and values of 29, Z59, and z, previously 
calculated into equation (19) and solve for the e,, . A common design period for 


residential and light commercial structures is r = 20 years (5 percent risk). 
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Fig. 16. (Lytton 1994) Edge Moisture Variation Distance for Center Lift Conditions (50 yrs) 
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3.3 EOGE LIFT CONDITIONS 
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Fig. 17. (Lytton 1994) Edge Moisture Variation Distance for Edge Lift Conditions (50 yrs) 
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Fig. 18. (Wray 1978) Edge Moisture Variation Distances for Edge/Center Lift Conditions (10yrs) 
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Equilibrium Suction 


Before the equilibrium suction profiles, h,,(z,), wet and dry limit suction profiles, 
hvin(Z,) and h,,,.(Z,), transient state suction profiles, h(z,,t), and depth to constant suction, 
z,, can be generated, the equilibrium suction , h,, must be known for particular soil 
profiles and surface conditions. Gay (1994) did extensive work in the area of 
calculating the mean volumetric moisture content for a given soil mass dependent upon 
the soil’s depth of available moisture, d,,,, and the location’s potential 


evapotranspiration. The depth of available moisture, d,,,, is defined as the maximum 


depth of moisture available for use by transpiring vegetation, which is stored within the 
soil zone down to the depth to constant suction. The Thornthwaite Moisture Index is 
used to ascertain if a climate is arid or humid (Thornthwaite 1948). Gay’s theory 
considers d,,, and TMI. With these, he developed a set of functional relationships that 
are used to calculate the mean volumetric moisture content, which are then used to 
calculate the equilibrium suction value for a particular soil profile and location. The 
approach he used to solve for these equations was developed by Juarez-Badillo (1975) 
in which functional domains for the problem are first established as shown in Fig. 19. 
Function (21) was developed to satisfy the boundary conditions in Fig. 19. Then 


through an assumption of linear proportionality of the rates of change of the two 


complete functions, Gay obtained, (22) and (23). 





1 | 
fay -(; =-¢| (21) 
aT df(d) 
YR = F(d) (22) 
a eee 
, re ae Se (23) 
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TMI + 60 =T 
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Fig. 19. (Gay 1994) Domains of Functions f(d) and g(T) 


Solving (23) for d gives (24). 


d = —=— (24) 


All symbols are defined in the list of symbols. Gay then determined how the 
coefficients y, d,, and T, depend upon d,,, using the pattern search technique and arrived 


at the expressions (25), (26), and (27). 


y = 0.0393374,,, + 1.357033 (25) 
d, = 0.449079d,,, + 0.304560 (26) 
T, = 0.062651d,,, + 59.53593 (27) 


These relationships facilitate the calculation of mean moisture depths, d,,, for all values 


of TMI for any value of available moisture depth between 10cm and 50cm, as illustrated 





a2 


in Fig. 20 (Gay 1994). The development of (24) is important because it is used to solve 
for the mean volumetric water content, 0,, The value 9,, is applied as a target to solve 
for h,, in an iterative process using Nieber’s (1981) equation (34). Nieber’s equation 
relates suction to volumetric water content. The method to apply this iterative process 


to solve for the equilibrium suction value can be found in Appendix A-90. 


Mean annual moisture depth (cms) 
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Fig. 20. (Gay 1994) Theoretical Relationships for Mean Annual Moisture Depth vs TMI 


The equation used to solve for the mean annual moisture depth, (24), developed 
by Gay is rearranged to arrive at the mean volumetric water content, (32). Below is the 
step by step process to obtain the equation to solve for 6,,. First note the following 
equations from Gay (1994), (28), (29), (30), (31). Fig 21 provides a graphical 
explanation for some of the symbols in these equations. All symbols are defined in the 


list of symbols. 


a (28) 
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im dry 
no Za (29) 
di 
a ce (30) 
a 
O ary at a (31) 
Dividing through equation (24) by z, gives (32). 
dam 
le 32 
= (32) 
ae - d, 
laa ae, 
ay 
| T, 
Solving for d,/z, in (29) and substituting in (32) gives (33). 
Ou = D ary 
0, = =F Oe, (33) 
da d, 
1+ 


The steps to arrive at equations which solve for 0,,, 0,,,, and 9,, using Nieber’s 
equation (34) are given below. 


] 
hh, | = (2-2. ] : 
mye e-8, 


(34) 
Introduce dummy variable r, (35). 








(35) 





34 


Solving (34) for 0 gives (36). 








aa 6.+r-0. 36 
een eS, 
B B 2 
: é Ih, | h,,,| hy) : ; 
Substituting a a or A back into (36) for r gives (37), (38), or (39) 


depending on which volumetric moisture content is to be found. 


MOISTURE CHANGES IN UNSATURATED SOILS 
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Fig. 21. (Lytton 1994) Moisture Changes in Unsaturated Soils Between Typical Wet and Dry States 
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5 ue (37) 


jens s A r (38) 





Me (39) 


The equilibrium suction value, h,, can then be solved through an iterative process using 
equations (33), (37), (38) and (39) as described in Appendix A-90. The equilibrium 
suction profile, h,,(z,) for any pervious soil profile is defined by decreasing the 
equilibrium suction value at the surface, h,, a centimeter in suction for every centimeter 
increase in depth (McQueen and Miller 1968). The process for determining the 


equilibrium suction profile is detailed in Appendix A-96. 


Suction Profiles 


It is important to note that moisture changes within a natural soil are determined 
by the ratio of the period of evaporation to the period of precipitation. A measure of 
this period of evaporation to the period of precipitation is determined through a method 
developed by Thornthwaite (1948). In the previous section the equilibrium suction is 
determined based on relations to the climatic water balance expressed in terms of the 
Thornthwaite Moisture Index (1948). When a marked separation occurs between wet 
and dry seasons, a large seasonal variation in soil moisture content occurs, whereas in 
areas which are either predominantly dry or predominantly wet for most part of the year, 
the changes in soil moisture content is not so marked. A location where predominant 
clays are abundant in which the climate is marked by this large separation in the periods 


of evaporation and periods of precipitation will experience problems with differential 
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movement to the greatest degree. The state of moisture is dependent on the state of 
suction and tends towards an equilibrium suction profile, as defined in the previous 
section, with a moisture source either at the boundaries or within the soil mass. 
When the state of suction is measured at intervals of depth in a soil profile, the resulting 
relation between suction and depth is termed the suction profile. Mitchell measured 
actual suction profiles in the area of Adelaide, South Australia and found a very 
consistent and typical trend that could be modeled. In his study he developed a method 
based on constitutive soil properties to model the suction profile within a soil profile 
base on the local annual weather cycle, n, soil diffusion, «(z,), equilibrium suction 
value, h,,, and the amplitude suction value, U,, The maximum limit suction profile, 
Hnax(Z;), and minimum limit suction profile, h,,;,(z;) begin at the limiting field condition 
suction value at the surface, and decrease at an exponential rate toward the equilibrium 
profile, h,,(z,). These limit suction profiles will approach the equilibrium suction profile 
at the depth to constant suction, z,,, which will be defined in the next section. The limit 
suction value depends on the type of surface conditions and is defined in the section 
titled “Typical Suction Levels in the Field”. The dry limit suction value is usually 4.5 
pF for soil with vegetation at the surface and 6.0 pF for bare soil at the surface. Refer to 
Fig. 22 for a typical maximum limit suction profile. The wet limit suction value, which 
is somewhat drier than the field capacity suction value, is usually 2.5 pF. Refer to Fig. 
23 for a typical minimum limit suction profile. Mitchell uses an idealized model of the 
annual weather surface cycle which he represents by (40). This is a good representation 
for design because it assumes that soil at the surface approaches the maximum and 
minimum moisture conditions at two distinct and separate times of the year. As stated 
earlier, this separation between wet and dry periods produces worst case differential soil 
movement scenarios. Mitchell combines the soil suction relationship with the 
relationship that describes the local surface weather cycle and arrives at equation (42). 
Mitchell found that the suction change due to the effects of climate, drainage and 


site cover is a periodic function of time. The suction at any time and at any depth in the 
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soil profile is determined by solving the diffusion equation, (41), for a particular 


Suction (pF) 


ot = 


Ul ui 
o © 





Fig. 22. Typical Dry Limit Suction Profile for Bare Soil at the Surface 


Suction (pF) 


0 

35 

70 

105 

gS 140 


3 175 


> 210 


J 245 
280 
$15 


350 4)" Equilitriun 
3385 Wet Linit 


Fig. 23. Typical Wet Limit Suction Profile for Bare Soil at the Surface 
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set of boundary conditions. The solution to arrive at the soil suction profile equation 


(42) can be found in Mitchell’s 1981 report (values of h(z,,t), h,,(z;,t), and U, are all in 


units of pF). 
A(z =.0,t) = cost ramet (z= Olu. 40 
oh oh 
Ot =O az" ( ) 





“Cex n-t 
NC Ze tien (7 yieluc co 2e-0-t (2) ae | (42) 


Based on the surface and depth boundary conditions, referred to as moisture effect 
cases, (42) is applied to produce suction profiles for many surface conditions and soil 
profiles. Appendix A contains detailed algorithms that describe precisely how to apply 
(42) to arrive at suction profiles resulting from combinations of different moisture 
effect cases and soil profiles. Later sections in “The Procedures, in a Software Design 
Format” chapter will provide a more detailed discussion with respect to application of 


(42) in producing suction profiles. 


Depth to Constant Suction 


The depth to constant suction is the depth in a soil profile to which there is no 
longer a significant seasonal change in soil suction. This depth has been reported to be 
anywhere from 6 feet to as much as 26 feet dependent upon the type of soil, seasonal 
weather trends, and the type of surface conditions such as vegetation and slope. Mohan 
and Rao reported where the active zone extended to a depth of 13 to 16 feet, in the 
black cotton soils of India. Kassiff reported movement of an exposed clay at depths up 
to 20 feet in Israel, with the largest differential movements being measured within an 


active zone of 6 % feet, where moisture changes are maximal (Mitchell 1980). As one 
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uses the procedures contained herein to predict the differential movement in soils, it 
becomes clear that most of the movement does in fact occur nearer the surface. This 
reduction in the amount of differential soil movement with an increase in depth is 
attributable to the combination of small seasonal changes in soil moisture content with 
an increase in depth, and the increase in restriction to volume strain due to overburden 
stresses with an increase in depth. The algorithm to solve for the depth to constant 
suction is based primarily on Mitchell’s equation (42) that is described in the previous 
section. 

Lytton and the author are responsible for developing the following steps to 
atrive at the equation to calculate the depth to constant suction. Only the limit profile 
aspect of Mitchell’s equation, (42), will be used. Setting the cosine term equal to one 
forces this equation to generate only the maximum and or minimum suction envelopes. 


Refer to (43). a 


(43) 


Mitchell’s equation in the form of (43) is used in most of the algorithms to calculate the 


—(Z;) ALS 
h(z,,t)=h,(z,)+Uje 8% 





volume strain for post equilibrium moisture effect cases. Equation (43) is rearranged to 


isolate the exponent term, (44). 





Zi 


= x 
h(z,,t)— h_(z.) = U-e ote (44) 
(44) is then divided through by U,, (45). 


nn 


a (Z;) 





h(z;,t)—h,,(Z) _ 
‘i = 


O 


-(z;) 
=e 


(45) 
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The last step is to take the natural log of both sides and solve for z,, (46). 





- nto 1 a 
U, n-t 
a(Z; ) 


The depth to constant suction, z,,, has been obtained when the amplitude at depth z, is 
equal to 0.2 pF. Later sections in “The Procedures, in a Software Design Format” 
chapter will provide a more detailed discussion with respect to application of (45) and 


(46) in calculating the depth to constant suction. 


Horizontal Velocity in Expansive Soils 


The ability to calculate horizontal velocity of water flow in soils is critical for 
the application of the soil volume strain theory. Predicting vertical differential soil 
movement when a horizontal barrier 1s employed, or when a moisture effect case is 
located a distinct distance from the slab’s edge depends upon the soil’s horizontal water 
flow characteristics. The horizontal velocity is calculated and then used to calculate a 
suction profile a distinct distance from another known suction profile. Lytton and the 
author developed the horizontal velocity equation by applying Darcy’s law, (47), to the 
unsaturated permeability relationship for clay soils developed by Mitchell, (18), to 
arrive at the horizontal velocity equation for unsaturated soils, (52). The steps to derive 


(52) are described below. First Darcy’s law, (47), is solved for dh, (48). 


dh 
Weay=-e( on 
dh = Ee dx (48) 


Darcy’s permeability constant is set equal to Mitchell’s unsaturated permeability, (49). 





4\ 


pep oadlecd 49 
= ; (49) 
Substitute (49) in for Darcy’s permeability in (48) and separate variables dependent on 
distance and variables dependent on suction to arrive at (50). 
Ne (Z, ) d bot dh 50 
ie (Z;) " h, i h 
Integrate the left side of (50) from zero to a distance x, and the right side of (50) from an 


entry suction value, h,,,,,(Z;), to an exit suction value, h,,;(z,). Refer to equation (51). 


V Ae) fax - Nexit (24) ib 


Sh 
_K o (2; i gh 0 Mentry (Zi yh 
Solve equation (51) for V,(z,) to arrive at equation (52). 
K,(Z,)h Sear) 
Z,) = a >2 
( hee (Z;) : 


The entry suction value is the suction value at x = 0. The exit suction value is the 
suction value at x. Once the horizontal velocity is known between two distinct locations 
in an incremental soil layer, the suction value for any point between those two distinct 
locations is solved for using equation (53). Equation (53) is derived by solving for the 


exit suction in (52). 


- V2; )X 
K, (Zz; )-he 
exit (Z; ) = ae (Z;) e (53) 


Post Construction Theory 


The worst case differential soil movement may occur in the first dry season 
subsequent to the construction of the slab during the wettest time of the year, or 
during the first wet season subsequent to construction of the slab during the driest time 
of the year. This type of moisture effect case is referred to as the post construction 


case. The post construction case is a transient case for calculating volume strain of 
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expansive soils immediately after construction and before the soil under the center of 
the slab has reached an equilibrium moisture content. For example, when a slab is 
constructed the soil is “capped” and the moisture content changes slowly toward its 
equilibrium value while the perimeter soil shrinks and swells with the change in 
seasonal moisture. This type of scenario provides conditions for the post construction 
moisture effect case. The moisture content of the “capped” soil looses or gains 
moisture for some period of time until the soil under the center of the slab reaches an 
equilibrium moisture content. The rate at which the soil approaches equilibrium is 
dictated by the soil’s diffusivity and unsaturated permeability. Once the soil under the 
slab reaches the equilibrium moisture content, subsequent suction distributions and 


differential movements must be determined using the post equilibrium methodology. 


The “Damped” Slab Suction Profile 


Lytton and the author developed the following set of equations to model the soil 
under the slab subsequent to construction in a transient case in which a soil suction 
profile approaches an equilibrium suction profile. The idea is to calculate the soil 
suction profile under the slab at the time of construction, hy,4(Zjsteonsvary)> and then 
dampen this profile until it approaches the equilibrium profile. The soil suction profiles 
under the slab and at the slab edge are solved, with respect to time, based on some 
knowledge of the number of annual weather cycles. The annual weather cycle will be 
assumed to vary according to the function described by (40). To demonstrate this 
method, the annual weather cycle is assumed to be one, n= 1. This means that, in one 
year there will be one wet season and one dry season. The driest, or wettest, month of 
the year must be known. In this case, June is selected as the driest month of the year. 
Fig. 24 is an example of a weather cycle modeled using (40) for bare soil at the surface 
(n= 1, and June is the driest month of the year). 

The approach used to derive the equations to damp the suction profile at time of 


construction back to the equilibrium suction profile was developed by Juarez-Badillo 
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(1975) in which functional domains for the problem are first established as shown in 


Fig. 25. Function (54) was developed to satisfy the boundary conditions in Fig. 25. 


For Bare Soil At The Surface 


Suction @ Surface (pF 
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Fig. 24. A Weather Cycle Modeled by a Cosine Function, n = 1, June is the Driest Month 


l | 
fih)= ay an (Case 1, when h,, is greater than h,,,,..) (54) 


m const 





] l 
f(h) = h. —h a (Case 2, when h,, is less than h,,,,,..) (54) 


m const 





The suction value with respect to time and depth for a profile that is modeled to tend 
toward an equilibrium suction was solved for cases 1 and 2. Both cases arrived at the 
same expression (59). The steps to arrive at (59) are described below. Through an 
assumption of linear proportionality of the rates of change of the function equations 


(55) is established, and the constant of proportionality, B, is equal to 1. 


df(h) dt 
£(h) roe (55) 
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increases 





Fig. 25. Domains of Functions h, t, and f(h) 


The left side of equation (55) is integrated from an arbitrarily chosen suction value, h, to 


a suction value half way between h,, and h,,,,, h, The right side of equation (55) is 


const? 


integrated from a time corresponding to h,, t,, to the time required to approach h, t. This 


integration is given by equation (56). 


= t 
m const B+] ; @ 6) 
hee— hye h. —h 


m const 


Equation (56) is rearranged to collect all constant variables on one side of the equation 


and set equal to a variable, r. Refer to equation (57). 


[B= Ho -(2 te fj - (57 
h -h) \h -hJ (tJ 


Now equation (57) is rearranged to solve for the suction value, h, given by (58). 


h _ ere cteniealn 


l+r (28) 
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Substitute (57) back in to (58) to arrive at the equation used to solve for the suction 
profile under the slab some time, t, after the slab was “capped” as a result of slab 
construction, (59). 


h, Sales -f B 
eee (Zebeee iain (Z) h ey oe mt 


m l 


hip olh 
1+] ———"]-t, ®t? 
at : 


h(z, t) = (59) 


The Time Equation 


Before equation (59) can be applied, a relationship must be derived to determine 
the amount of time required for water to move upward or downward to change the 
suction at any given depth half way between its initial value and the equilibrium value 
of suction. This time equation for t, is given by (61). The derivation of (61) is as 
follows. Darcy’s law of permeability is the relationship used to introduce the time 
variable. Velocity is defined as a distance through which water flows in soil, x, traveled 
in some amount of time, t,. Velocity described in terms of permeability of water in 
expansive soils in which the time, t,, required to obtain half of the change in suction (the 


time to obtain half of the change in suction was arbitrarily chosen), 0.5- Ah, is given by 


(60) 
Nt 7 - pe.) ? = (60) 


D4 





Equation (60) is solved for t,, (61) In equation (60), the variable x is replaced with 


equation (44), where the ratio of suction values is set equal to one half. 


i gale : (05) 
Nn: 7 


eimai Si le (61) 
09 9.4343 -(h(Z; st gonst) — Bin (Z)) 


The time variable represented by (61) is used in (59) to solve for the suction profile 


t, = 


some time, t, after soil is “capped” as a result of slab construction. 
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THE PROCEDURES, IN A SOFTWARE DESIGN FORMAT 


Introduction to the Software Design Approach 


The tools described in “The Theory of Differential Soil Movement” chapter are 
all that is needed to develop the equations and terms used in the procedures described in 
this chapter. Appendix B is a graphic representation suggested as a flow of vertical 
differential soil movement algorithms in a windows environment type software design 
layout. The suggested software design layout, Appendix B, is modeled after the 
PTISLAB WIN 1.0 software package owned by Geostructural Tool Kit, Inc. (1996). 
The focus of this document is not software design or marketing, those skills are left to 
experts in that area of work. The focus 1s to present the equations and terms in easy to 
follow algorithms, as described in Appendix A, so that engineers or programmers can 
apply this soil volume strain technology in everyday practice. 

A question that a practicing engineer may ask 1s, “What practical engineering 
cases can these procedures address?” The answer is, cases involving a surface drained 
soil profile in which the slope is stable and surface conditions are those described 


below: 


1. Bare soil at the surface. 
Grass at the surface with shallow roots. 


A flower bed at the surface with a known depth of flower bed zone, d 


zone’ 


ae 


A tree with a known depth of root zone, d,,,..- 
The procedures consider design effects that are introduced to change the magnitude of 
the vertical differential soil movement. These design effects include: vertical barriers, 
horizontal barriers, and the effects of locating the surface moisture effect a distinct 
distance from the foundation, d,,,,... All of these design effects can be rigorously 


analyzed to enable an engineer the ability to choose the effect which will provide the 
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desired design results. 

The algorithms contained herein address all cases and design effects listed above 
by way of two design approaches: (1) the post construction approach, and (2) the post 
equilibrium approach. An engineer would apply the algorithms, used in both 
approaches, to address all possible combinations of cases and design effects presented. 
The engineer compares the results of these cases to find scenarios that produce vertical 
differential soil movement that generates maximum bending, shear and deflection. This 
information is used to generate adequate foundation designs. The most promising 
scenario is to use the power of today’s computing technology to perform the iterative 
calculations and tests required to arrive at these design quantities. A programmer could 
apply software technology, such as PTISLAB WIN 1.0, to model a slab’s interaction 
with the differential soil movements produced by the algorithms contained herein to 
arrive at these maximum design cases. 

The algorithms are presented in a format that can be easily programmed using 
simple spreadsheet type software. First, input required for all calculations are entered. 
Soil properties and other related soil characteristics are calculated and stored for use in 
algorithms presented later. Equilibrium suction values, h,, and depth to constant 
suction values, z,,, are calculated for a specific location and specific sets of surface 
conditions. With the above complete, all moisture effect cases and design effects listed 
previously are addressed using the post construction solutions and the post equilibrium 


solutions. 


Input 

When considering all cases and all calculations required to complete algorithms, 
an accounting of the variables required for input becomes important. Great effort is 
made to keep the complex nature of data required for input down to a minimum. The 
goal, relative to input of data, is to ensure that all variables required for input are items 
that are easy and economical for practicing engineers to obtain. This is the same goal 


authors had in mind while developing the simplified method for identifying 
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predominant clay minerals (Pearring 1968, Holt 1969 and Mojekwu 1979), the methods 


for obtaining the suction compression index and the mean principal stress compression 


index (McKeen 1981 and Viayvergiya and Ghazzaly 1973), and the method for 


obtaining the edge moisture variation distance (Jayatilaka et al.1992). Below is a full 


list of items required for input to perform every calculation in the algorithms described 


in Appendix A. 


A,B 
B 


Ie 
M, 


S 


PL, 
TMI 
Xp 
Yb 
d 


am 
dibar 
Gls sires 
as 
zone 
c 


Ray 


Gardner’s coefficients (Appendix A-90 and Table 2) 

Number of boring samples (Appendix A-84) 

Depth to the c™ suction measurement in a sample (Appendix A-87) 
Depth to the bottom of layer s (Appendix A-84) 

Liquid limit for layer s, in percent (Appendix A-84) 

Number of layers in the b" boring sample (Appendix A-84) 

Number of suction measurements in a sample (Appendix A-87) 
Plastic limit for layer s, in percent (Appendix A-84) 

Local annual Thornthwaite Moisture Index, (Fig. 26 and Appendix A-83) 
X coordinate for the b" boring sample (Appendix A-84) 

Y coordinate for the b" boring sample (Appendix A-84) 

Available moisture depth of the soil mass (Appendix A-90) 

Size of horizontal barrier (Appendix A-86) 

Distance effect is away from the edge of barrier/slab (Appendix A-86) 
Depth of vertical barrier (Appendix A-86) 

Depth of the root or flower bed zone (Appendix A-86) 

c suction measurement in a sample (Appendix A-87) 

Dry suction limit for soils in the field, bare soil at the surface usually 
has a value of 6.0pF and soil with vegetation at the surface usually 
has a value of 4.5pF (It may be requested in various algorithms. 


Appendix A-90 is the first time hg,, 1s requested as input). 
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hy Wet limit suction for soils in the field, usually 2.5pF (It may be requested 
in various algorithms. Appendix A-90 is the first time h,, is used). 
n Local annual weather cycle (Appendix A-83). 


t The month the slab was constructed (Appendix A-107) 


const 


te Driest month of the year (Appendix A-107) 

f. The month the soil volume change is sought (Appendix A-107) 
LL. Percent passing 2 micron for layer s (Appendix A-84) 

V. Percent passing number 200 sieve for layer s (Appendix A-84) 


More detailed explanations of how these variables are used 1s described in the following 
sections of this document. Appendix A gives a complete process flow detailing the use 
of the variables listed above and the variables listed in the list of symbols. 

Input, as requested by the process flow in Appendix A and as graphically 
represented in the software design layout sheets in Appendix B, is broken down into six 
major categories: 

1. Project Information input. This includes project name, date the project started, 
name of the engineer assigned to the project, and the site location. This information 
is basic information an engineering firm may want to establish in a database so that 
information can be readily accessed. The only information critical to calculating 
differential soil movement is the site location. The site location can be used as the 
means to locate the TMI for a site without requiring the user to manually research 
this bit of information. Refer to Appendix A-82 for a process flow and Appendix B- 
117 for a program design layout. The TMI for a particular site is found in Fig.26. 

2. Foundation Layout. This is basic surveying information to establish a three 
dimensional grid for logging information related to appropriate physical locations. 
Refer to Appendix A-83 for process flow and Appendix B-118 for program design 
layout. 

3. Soil Layer Property Input. These values include boring sample data and soil layer 


property data. The data entered here will provide enough information to generate a 
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three dimensional picture of the soil profile. Refer to Appendix A-84 for a process 
flow and Appendix B-119, 120, 126 for program design layouts. The soil profile 


used in all example calculations is shown in Table 1. 
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Fig. 26. (Wray 1978) Thornthwaite Index for Texas (20 year average 1955-1974) 


Table 1. Soil Profile Considered in all Sample Calculations 


Soil Va 
Layer No. Ibs/ft° 
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4. Moisture Effect Input. This gives the user the option to input the type and location 
of the moisture effect cases considered in vertical differential soil movement 
calculations. Refer to the plan view, Appendix C-132, for a view of the foundation 
used as the example throughout this document. The input for location and moisture 
effect case type can be found listed on the program design layout Appendix B-121. 
For example, case “1-5” describes the sidewalk with grass covered lawn on the 
north and west portions of the building shown in Appendix C-132. The number “1”, 
in “1-5”, describes this case as the first sequential case of this type. The number 
“5” in “1-5”, describes this case to be grass at the surface with a horizontal barrier. 
Additional information 1s shown in Appendix B-121 and is described in the process 
flow sheets Appendix A-85, 86. 

5. Measured Suction Profile Input. Geotechnical engineers may identify some 
unique moisture effect cases that cannot be modeled in the procedures to address the 
typical moisture effect cases contained herein. Those special cases are addressed 
using initial and final measured suction profiles. The input of these measured 
suction profiles is described in the process flow sheet Appendix A-87, and shown in 
the program design layout sheet, Appendix B-122, 123. 

6. Time Input for Post Construction Cases. Algorithms for post construction 
solutions are time dependent. Calculating the post construction differential soil 
movements can only be accomplished if the driest (or wettest) time of year is 


known, t,,,, the month of year the foundation was built is known, t and the month 


const? 


the soil volume change is sought is known, t,.. Input of these quantities is described 


in the process flow sheet Appendix A-107. 


Calculated Soil Properties and Soil Characteristics 
Once the liquid limit, the plastic limit, the percent passing the number 200 sieve, 
the percent smaller than two microns, and the dry unit weight is entered for the soil 


profile considered, the remainder of the required soil properties and characteristics are 
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solved. The process flow sheet, Appendix A-88, describes the algorithm by which these 
remaining soil characteristics are solved. Definitions for all symbols are found in the 
list of symbols. 

The plasticity index, PI(z,), 1s the difference between the liquid limit, LL(z,), and 
the plastic limit, PL(z,). The percent fine clay, Clay(z,), is calculated by dividing the 
percent smaller than two microns, 1(z,), by the percent passing the number 200 sieve, 
v(z,). The activity ratio, Ac(z,), 1s calculated by dividing the PI(z,) by Clay(z,). The 
cation exchange capacity, CEC(z,), is the larger of the values PL(z,)'""’ or LL(z,)°?”. 

The cation exchange activity, CEAc(z,), is calculated by dividing CEC(z,) by Clay(z,). 
Using Fig. 15, the matrix suction compression index, ¥;,;99, 18 the value found at the 
intersection of CEAc(z,) and Ac(z,). To obtain the y,(z,) simply multiply Clay(zi) by 
Ynioo: Lhe slope of the suction versus the gravimetric water content line, S(z,), is 


approximated by (16). Mitchell’s diffusion coefficient, «(z,), is approximated by (17). 








Mitchell’s unsaturated permeability, K,(z,)-|h,], 1s calculated by (18). Using Fig. 16 
for center lift conditions and Fig. 17 for edge lift conditions, the edge moisture variation 
is found by finding the curve that corresponds with the soil's diffusion coefficient, a(z,), 
then finding the intersection made by drawing a vertical line corresponding to the local 
TMI (for locations 1n Texas the local TMI can be found in Fig. 26), and finally drawing 
a horizontal line from that intersection across to the y axis identifies the appropriate 
edge moisture variation distance. For each case of center lift and edge lift, the largest 
edge moisture variation value, considering all incremental soil layers, is the soil 


profile’s edge moisture variation distance, e,,, and e,,.. 


Procedures to Calculate the Equilibrium Suction 


The next step in the process of calculating the vertical differential soil movement 
is to establish a baseline equilibrium suction, h,, (units of pF). All suction profiles are 
calculated relative to this baseline. The equilibrium suction, h,, 1s calculated for the soil 


layer located at the depth to constant suction. The depth to constant suction, z,,, is 
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dependent upon the h,, Therefore, if the depth to constant suction, z,,, does not fall 
within the soil layer in which the h,, was calculated, then the h,, must be re-calculated 
for the soil layer located at the depth to constant suction. This is an iterative process. 
The process flow sheet which describes the procedure for calculating h,, is found on 
Appendix A-90. The procedures for calculating z,, are found in the next section. The 
equilibrium suction value is dependent upon the following inputs: TMI, d,,,, hy, hyy, 9,, 
0, A, and B. For the state of Texas the TMI is obtained from Fig. 26. Values for 0, 9,, 
A, and B for a particular classification of soil are listed in Table 2. As reported by Gay 
(1994), a conservative estimate for the d,,, can be taken as 30 cm for all expansive soils 
in which the d,,, 1s not known. All variables have default values. Refer to Appendix A- 
90. Once all input is entered y, d,, T,, 0,,, 94,, are calculated using (25), (26), (27), (38), 
and (39) respectively. These values are used in equations (33) and (37) to solve for h,,. 


A value of h,, 1s picked and an iterative process for solving for 9,, is performed until 


equal values of 0,, satisfies both equations (33) and (37). 


Depth to Constant Suction 


The depth to constant suction is the depth in a soil profile to which there is no 
longer a significant seasonal suction change. The definition of no significant seasonal 
suction change is 0.2 pF from Lytton’s 1994 paper. Depth to constant suction is 
dependent upon the diffusivity of the soil, the number of local annual weather cycles, 
and the suction amplitude at the surface. As discussed in the “Soil Volume Strain” 
section, the soil’s ability to change volume is dependent upon the seasonal change of 
moisture and the magnitude of the mean principal stress present to suppress the volume 
change. Using Mitchell’s equation (42) there will always be a reduction in the suction 
amplitude with an increase of depth. Additionally, an increase in depth introduces 
added restraint to soil volume strain due to the effect of overburden. Considering these 
two aspects of the soil volume change equation it becomes apparent that for depths 


greater than 30 feet there is usually no longer any significant changes in suction. The 





depth to constant suction, z,,, will be used as the depth beyond which there is no 


significant change in suction. 


GM-GC1 GM-GC 


GM1 GM 
GM2 GM 
GM3 GM 
GM4 GM 
GP 1 GP 
GW1 GW 
GW2 GW 
GW3 GW 
SM-SC1 SM-SC 
SM1 SM 
SM2 SM 
SM3 SM 
SM4 SM 
SM5 SM 
SM6 SM 
SM7 SM 
SM8 SM 
SM9 SM 
SM10 SM 
SM11 SM 
SM12 SM 
SP-SM1 SP-SM 
SP-SM2 SP-SM 
SP1 SP 
SP2 SP 
SP3 SP 
SP4 oP 
SRS SP 
SW-SP1 SW-SP 
Swi SW 
SW3 SW 
lr CL 
CL2 CL 
CL3 CL 
CL4 CL 
CL5 CL 
CL6 CL 
CL7 CL 
CL8 CL 


Unified 
Soil 


Table 2. (Lytton et al. 1989) Gardner’s Coefficients 


0.42 
0.259 
0.382 
0.451 
0.608 
0:2595 
0.443 
0.506 
0.252 
0.656 
0.396 
0.619 
VAg2Z 
0.655 
0.623 
1.024 

OFZ 
0.831 
0.721 
0.818 
1.132 

2.16 
0.574 
0.583 
0.712 
0.736 
0.818 
0.593 
0.583 
0.773 
0.551 
0.385 
1.183 
0.621 
0.785 

0.61 
0.825 
0.964 
0.808 
Org 


0.004 
0.152 
0.04 
0.066 
0.043 
0.065 
0.039 
0.596 
0.309 
0.013 
0.016 
0.001 
0.11 
0.039 
0.015 
0.011 
0.023 
0.208 
0.01 
0.018 
0.029 
0.042 
0.048 
0.095 
0.0554 
0.024 
0.042 
0.076 
0.053 
0.012 
0.162 
0.111 
0.065 


0.004 


0.018 


0.028 


0.637 
0.269 
0.648 
0.251 
0.478 
0.56 
0.302 
0.318 
0.319 
Oia 
0.562 
1.023 
0.339 
0.468 
0.835 
0.671 
0.549 
0.436 
0.835 
0.806 
0.745 
0.501 
0.769 
0.613 
0.79 
0.951 
0.9 
0.665 
0.809 
1.082 
0.5852 
0.616 
0.417 
0.976 
0.61 
0.957 
QisZa 
leon 
1.593 
0.344 


Porosity 


0.296 
0.206 
0.276 
0.311 
0.378 
0.203 
0.307 
0.416 
0.355 
0.396 
0.284 
0.382 
0.544 
0.396 
0.384 
0.506 
0.419 
0.454 
0.419 
0.45 
0.531 
0.684 
0.365 
0.368 
0.416 
0.424 
0.45 
0.372 
0.368 
0.296 
0.355 
0.278 
0.542 
0.383 
0.44 
0.379 
0.452 
0.491 
0.447 
0.493 


0.0296 
0.0206 
0.0276 
O10S 10 
0.0378 
0.0203 
0.0307 
0.0416 
0.0355 
0.0396 
0.0284 
0.0382 
0.0544 
0.0396 
0.0384 
0.0506 
0.0419 
0.0454 
0.0419 
0.045 

0.0531 
0.0684 
0.0365 
0.0368 
0.0416 
0.0424 
0.045 

0.0372 
0.0368 
0.0296 
0.0355 
0.0278 
0.0542 
0.0383 
0.044 

0.0379 
0.0452 
0.0491 
0.0447 
0.0493 


1083.777 
12.2224 
111.1578 
2/.00574 
69.90875 
94.58668 
51.39672 
3.489424 
6.745925 
452.9567 
227.9712 
10543.87 
19.843 
75.32435 
455.9411 
426.194 
153.9119 
13.12009 
683.9116 
355.4082 
191.6911 
75.46589 
122.3936 
43.17938 
111.2987 
372.2106 
189.1258 
60.83961 
121.5414 
1006.512 
Ao TONG 
37.21149 
40.1871 
0 
1018.451 
0 
185.2369 
0 
0 
78.85731 





0.637 
0.269 
0.648 
0.251 
0.478 
0.56 
0.302 
0.318 
0.319 
On 
0.562 
1.023 
Oe2 
0.468 
0:639 
0.671 
0.549 
0.436 
0.835 
0.806 
0.745 
0.501 
0.769 
0.613 
0.79 
0.951 
0.9 
0.665 
0.809 
1.082 
0.5852 
0.616 
0.417 
0.976 
0.61 
0.957 
0.523 
1.91 
1.593 
0.344 
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CL9 CL 
CL10 CL 
CL11 CL 
ML8& ML 
ML9 ML 
ML10 ML 
ML11 ML 
ML12 ML 
lE{2 CL 
ML-OL1 ML-OL 
ML-OL2 ML-OL 
ML-CL1 ML-CL 
ML-CL2 ML-CL 
ML-CL3  ML-CL 
ML1 ML 
ML2 ML 
ML3 ML 
ML4 ML 
MLS ML 
ML6 ML 
ML7 ML 


The depth to constant suction varies dependent upon the soil strata, the local 


Table 2. Continued 


AWL XWL 

0.562 6.9E-06 1.551 
0.517 0.024 0.359 
0.415 0.064 0.247 
0.858 0.003 0.973 
0606 0.022 0.707 
2.163 0.038 0.61 
0.603 0.031 0.722 
0.821 eee oe 
0.972 0.004 0.474 
1.721 0.064 0.535 
0.885 O 1.362 
0.692 0.066 0.365 
0.642 0.013 0.634 
0.645 0.001 1.052 
1.015 0.003 1.054 
1.067 Deel 25i7, 
1.119 0.032 0.681 
1:558 401665". 0-411 

0.88 O 1.019 
0.825 O 1.474 
0783-50012 05/5 


Porosity 


0.36 
0.341 
0.293 
0.462 
0.377 
0.684 
0.376 
0.451 
0.493 
0.632 
0.469 
0.409 
0.391 
0.392 
0.504 
0.516 
0.528 
0.609 
0.468 
0.452 
0.439 


0.036 
0.0341 
0.0293 
0.0462 
0.0377 
0.0684 
0.0376 
0.0451 
0.0493 
0.0632 
0.0469 
0.0409 
0.0391 
0.0392 
0.0504 
0.0516 
0.0528 
0.0609 
0.0468 
0.0452 
0.0439 


9154077 
95.23328 
27.59434 
3132.411 

231.514 
107.2053 
170.0741 

0 
744.6291 
53.55747 
0 
35.11204 
331.1743 
111271.97 
3/774.668 
0 
149.9167 
39.6357 1 
0 
0 
313.1978 


1.551 
0.359 
0.247 
0.973 
0.707 

0.61 
0.722 
2.059 
0.474 
0.535 
1.362 
0.365 
0.634 
1.052 
1.054 
1.257 
0.681 
0.411 
UE 
1.474 
0575 


weather conditions, and the type of soil surface conditions. The depth to constant 


suction, z,,, will be calculated for the following cases of soil surface conditions, 


Appendix A-91: 


1. Bare soil at the surface (Appendix A-92). 


2. Grass at the surface with shallow roots (Appendix A-92). 


3. Flower bed at the surface with a known flower bed zone depth, d,,,,, (Appendix A- 


93). 


4. A tree at the surface with a known depth of root zone, d 


5. And measured suction profiles (Appendix A-94). 


zone 


(Appendix A-93). 
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The simplest case would be that of bare soil or grass at the surface with only one 


soil layer the entire depth of moisture active zone. For this case equation (46) is all that 
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is needed to solve for z,,. Simply replace the numerator in the natural log function with 
0.2 pF, substitute entered or previously calculated values for all other variables and 
solve for z,,. The suction amplitude, U,, is the difference between maximum and the 
minimum surface soil suction values. For bare soils the suction amplitude is U, = hy, - 
h,,, where h,,, = 6.0 pF. For soils with grass at the surface the suction amplitude is U, = 
gy - Hy, where hy, = 4.5 pF. Refer to the “Typical Suction Levels in the Field” for 
typical values of suction. 

In most soil profiles there is usually at least two soil layers to consider. In cases 
considering bare soil or grass at the surface, with multiple soil layers, the algorithms in 
Appendix A-92 are used to solve for z,,. The procedures are simple. 

1. First solve for the ratio of the suction difference between the top increment and 
bottom increment of a soil layer, P,, using equation (45). 

2. Multiply P, by the suction amplitude at the top of the layer which gives the suction 
amplitude at the bottom of the layer, U,.. 

Be vest U,.: 

a. If U,, 1s less than 0.2 pF, then the depth to constant suction is in this 
layer, in which case z,, 1S solved using (46). The initial U, is the same as 
the one described in the previous paragraph. All subsequent U,’s are 
solved using the procedure in paragraph 3.b. described below. 

b. If U,, is grater than 0.2 pF, then set U, equal to U,, and repeat these three 
steps. 

The depth to constant suction 1s the z,, calculated for the last layer added to the value of 
the depth to the top of this layer. 

The procedures for calculating the z,, for flower beds and trees 1s basically the 
same with two distinct differences: (1) The suction amplitude is the difference between 
the limit suction value and the equilibrium suction value. (2) The suction amplitude at 
the surface is constant to the depth of moisture effect zone, d,,,... For a flower bed at 
the surface the suction amplitude is the difference between the equilibrium suction value 


and the wet limit suction value, U, =h,, - h,,. For a tree at the surface the suction 
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amplitude is the difference between the dry limit suction value for vegetation and 
equilibrium suction value, U, = hy, - h,,. Refer to Appendix A-93 for a detailed process 
flow of these procedures. Refer to “Typical Suction Levels in the Field” and 
“Procedures to Calculate the Equilibrium Suction” for values of (h,, and hy,) and h,, 
respectively. 

When making calculations of the depth to constant suction, z,,, for the case of 
measured suction profiles, the depth to constant suction can be taken as the value of z,, 
calculated for one of the four previously explained cases if the surface conditions are the 
same. However, a check must be made to ensure that the z,, in a particular measured 
suction profile case does not extend beyond those values of z,, in the previously 
described cases. If the z,, extends beyond those values calculated using the procedures 
described in the previous four cases then the z,, must be calculated using the procedures 
described in Appendix A-94. The procedures are the same as those described above for 
bare soil and grass at the surface with one major difference. The difference is that the 
process 1s initiated at the bottom of the measured suction profile. This z,, value 
calculated for a measured suction profile is then compared to the z,, for the surface 
conditions most similar to the measured profile surface conditions. The depth to 


constant suction is taken as the larger of these values. 
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POST EQUILIBRIUM AND POST CONSTRUCTION SOLUTIONS 


At this point in the procedures all inputs and calculated values required to generate 
suction profiles for soil volume strain calculations have been explained. The procedures 
to generate soil suction profiles and the corresponding soil movement for moisture effect 
cases and design effects are explained using the following two procedural approaches: (1) 
the post equilibrium approach, and (2) the post construction approach. The post 
equilibrium approach refers to procedures to solve for differential movement of expansive 
soils after the soil under the slab has reached an equilibrium moisture content. When a 
slab is placed on an expansive soil, there is a time period in which the moisture content 
under the slab changes from the suction profile as it was when the slab was placed to 
approach the equilibrium suction value. The algorithms to solve for differential soil 


movement during this time period are handled by the post construction approach. 


Post Equilibrium Approach 


Once the soil under the slab reaches the equilibrium suction value, the differential 
soil movement is calculated using the differential in soil suction profiles between the 
equilibrium suction profile under the slab and the suction profile at the edge of the slab. 
The suction profile at the edge of the slab is dependent upon the seasonal moisture change. 
The restraint to soil volume increase at any depth is dictated by the amount of overburden 
and surcharge stress. The post equilibrium algorithms to solve for the limit cases are: (1) 
differential soil movement causing a maximum center lift distortion mode, and (2) 


differential soil movement causing a maximum edge lift distortion mode. 


Equilibrium Suction Profiles 


The first calculations that are made in the post equilibrium approach are the base 
line equilibrium suction profiles. The equilibrium suction profile is generated by reducing 
the equilibrium suction value calculated for a specific moisture effect case by one 
centimeter in suction for every centimeter of descent in soil depth. Care should be taken to 


convert the calculated equilibrium suction value from units of pF to units of centimeters. 
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Once the equilibrium suction profile is generated these values of suction in units of 
centimeters should be converted back to units of pF to make further calculations. Refer to 


Appendix A-96 for algorithms to calculate equilibrium suction profiles. 


Suction Profiles Causing The Center Lift Distortion Mode 


Center lift distortion modes are caused by a differential in the equilibrium suction 
profile under the slab and the dry limit suction profile at the edge of the slab. Once the 
equilibrium suction profiles are calculated for each case considered, the dry limit suction 
profiles for the slab edge moisture effect conditions are calculated. The dry limit suction 
profiles for bare soil and grass at the slab edge are calculated using Appendix A-97, using 
variables unique to the specific moisture effect case and location. The dry limit suction 
profiles for a flower bed and a tree at the edge of slab are calculated using the Appendix A- 
98. The suction profiles are calculated by the application of (43). The differences in these 
processes, dependent upon moisture effect case and location are: the surface suction values, 
the soil properties, the baseline equilibrium suction profiles, and depth to constant suction 
values. For all moisture effect cases, the suction amplitude change and the depth are 
initialized at every soil layer change. The cases of a tree and a flower bed at the surface 
have one other major difference, these cases include the depth of moisture effect zone. The 
suction profile within the tree root or flower bed zone is constant and equal to the limit 


suction value for the full depth of the moisture effect zone, d,,,.. Refer to Fig. 27 for an 


example of the dry limit post equilibrium suction profiles for a flower bed at the surface 
with a 4 foot moisture effect zone depth. Refer to Fig. 28 for an example of the dry limit 


post equilibrium suction profiles for a tree at the surface with 20 foot deep roots. 


Suction Profiles Causing The Edge Lift Distortion Mode 


Edge lift distortion modes are caused by a differential in the equilibrium suction 
profile under the slab and the wet limit suction profile at the edge of the slab. Once the 
equilibrium suction profiles are calculated for each case considered, the suction profiles for 
the slab edge moisture effect conditions are calculated. The wet limit suction profiles for 


bare soil, grass, and trees at the slab’s edge are calculated using Appendix A-99. The wet 
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limit suction profiles for a flower bed at the edge of the slab is calculated using the 


algorithm described in Appendix A-100. The differences in these processes which depend 
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Fig. 27. Post Equilibrium Suction Profiles for Flower Bed at the Surface (Dry Limit) 
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Fig. 28. Post Equilibrium Suction Profiles for a Tree at The Surface (Dry Limit) 
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upon the moisture effect case and location are: the soil properties, the baseline equilibrium 
suction profiles, and depth to constant suction values. Note that the suction profile is 
constant and equal to the limit suction value for the full depth of the flower bed zone, d.,,... 
The suction amplitude change and the depth are initialized at every soil layer change. 


Refer to Fig. 29 for an example of the wet limit post equilibrium suction profiles for a 


flower bed at the surface. 
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Fig. 29. Post Equilibrium Suction Profiles for a Flower Bed at The Surface (Wet Limit) 


Vertical Differential Soil Movement 

Expansive soil volume change 1s predicted by applying the theory explained in the 
‘Soil Volume Strain” section of the previous chapter. The moisture effect cases considered 
consist of: 
1. Bare soil at the surface, 


2. Grass at the surface with shallow roots, 





62 


3. A flower bed at the surface with a known depth of flower bed zone, d nd 


zone> & 
4. A tree at the surface with a known depth of root zone, d.,,.. 
The procedures also consider the application of design effects to provide the best possible 
conditions so foundation design performance is at a pinnacle. The design effects 
considered include: 

1. Vertical barriers, 

2. Horizontal barriers, and 

3. The effect of locating the moisture effect case a distinct distance from the slab edge. 
Explanation of the procedures to predict vertical differential soil movement is broken down 
into three categories: (1) all moisture effect cases with or without a vertical barrier, (2) 
moisture effect cases with or without a horizontal barrier, and (3) moisture effect cases 
located a distinct distance from the slab edge. 

Post equilibrium vertical differential soil movement can be determined for all 
moisture effect cases with or without a vertical barrier by simply applying the algorithms 
described by Appendix A-103. This process can be broken down into three basic steps: 

(1) All pertinent variables and information collected or calculated in previous procedural 
steps are retrieved for use in this process. (2) Tests are performed on variables to ensure 
the appropriate calculations and quantities are used, based on Lytton’s volume strain 
theory. (3) Incremental strain and corresponding vertical soil movement is summed up for 
every soil increment within the moisture active zone of the soil profile. The total vertical 
differential soil movement quantity is defined by the variable Az. Refer to Fig. 30 for an 
example of bare soil dry limit suction profiles (center lift distortion mode) for a slab edge 
that has a four foot deep vertical barrier. 

Post equilibrium vertical differential soil movement can be computed for all 
moisture effect cases, except the tree case, with or without a horizontal barrier by simply 
applying the algorithms described by Appendix A-104. The case of the tree is unique in 
that the tree roots will be active regardless of whether a horizontal barrier exists or not. If 
the tree is killed because of the placement of the horizontal barrier then the point becomes 


mute. The process in Appendix A-104 is broken down into four basic steps: (1) All 
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Fig. 30 Bare Soil Dry Limit Suction Profiles for a Slab with a 4ft Vertical Barrier 


pertinent variables and information collected or calculated in previous procedural steps 

are retrieved for use in this process. (2) The horizontal velocity for each incremental soil 
layer is calculated based on the edge moisture variation distance and the size of horizontal 
barrier. This step is based on equation (52). (3) The suction at the slab edge is calculated 
using equation (53). (4) This slab edge suction value, h,,,,(Z,), is used in the place of the 
final suction value in Appendix A-103 to make vertical differential soil movement 
calculations, Az. Refer to Fig. 31 for an example of wet limit suction profiles (edge lift 
distortion mode) for a four foot deep flower bed at the edge of a four foot horizontal barrier 
constructed adjacent to the edge of the slab. 

If a flower bed or tree is located a distinct distance from the slab edge, the vertical 
differential soil movement is calculated by applying Appendix A-104. These are the only 
two cases that can be located a distance from the slab and still have an effect on vertical 
differential soil movement under the slab. The differential movement caused by the tree or 


flower bed located a distance from the slab edge should be compared to the differential 
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Fig. 31. Suction Profiles for Flower Bed at the Slab Edge with a 4 ft Horizontal Barrier 


movement caused by the surface condition that exists between the slab edge and the 
tree/flower bed. The maximum differential soil movement produced by either the 
tree/flower bed case or the conditions between the slab edge and the tree/flower bed case 
controls. This process is broken down into four basic steps: (1) All pertinent variables and 
information collected or calculated in previous procedural steps are retrieved for use in this 
process. (2) The horizontal velocity for each incremental soil layer 1s calculated based on 
the edge moisture variation distance and the distance the moisture effect is away from the 
slab edge, d.jce. Lhis step is based on equation (52). (3) The suction at the slab edge is 
calculated using equation (53). (4) This slab edge suction value, h,,,.(Z,), is used in the 
place of the final suction value in Appendix A-103 to make vertical differential soil 
movement calculations, Az. Refer to Fig. 32 for a dry limit suction profiles (center lift 
distortion mode) for tree roots that are 20 ft deep and located four feet away from the edge 


of a slab with a four foot vertical barrier. 
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Fig. 32. Suction Profiles for a Tree 4ft From the Slab Edge with a 4ft Vertical Barrier 


Post Construction Approach 


As described 1n the “Post Construction Theory” section earlier, the worst case of 
vertical soil movement in expansive soils may occur in the post construction time period of 
the life of a structure. Lytton and Woodburn’s 1973 paper provides a perfect illustration of 
this worst case scenario occurring in the post construction time period during the life of a 
school building. The procedures to apply the post construction theory is basically a two 
step process. First, all suction profiles considered are calculated. Second, these suction 
profiles are used to calculate vertical soil movement at the center of the slab, Az,..,, and 


combined with the vertical soil movement at the slab edge, Az,,.,., to arrive at a total 


edge? 


vertical differential soil movement quantity, Az,,,.1- 


Post Construction Suction Profiles 


The suction profiles must be calculated for two time periods: soil suction profiles 


at the time of construction, t,,,,,, and soil suction profiles at the time the vertical differential 
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soil movement is sought, t,,.. Soil suction profiles at the edge of the slab are generated for 
the time at which the vertical differential soil movement is sought. Soil suction profiles 
beneath the slab are generated at the time of construction and at the time the vertical 
differential soil movement is sought. 

Appendix A-107 is the process flow sheet that: describes the time variables, lists a 
step by step process for entering the required times, and lists the steps to calculate the time 
periods needed to make soil suction profile calculations. Appendix B-127 is a picture of 
the program design layout sheet for the time input and display of suction profiles both 
under the slab and at the edge of the slab for time that transpires from the driest month to 
the month the volume change 1s sought. In Appendix B-127, the suction profile under the 
slab at the time of construction is the thin solid line. The dashed line 1s the suction profile 
under the slab six months after construction. Notice that this dashed line 1s tending toward 
equilibrium. The equilibrium profiles are the vertical thin lines with tick marks. Finally, 
the thick solid line is the suction at the edge of the slab with a vertical barrier. To generate 
suction profiles illustrated in Appendix B-127 the following time inputs are required: (1) 


and (3) the 


const? 


the driest month of the year, t,,, (2) the month the slab is constructed, t 
month the soil volume change is sought, t,,. These entered values of time are required for 
two primary purposes: (1) to establish a baseline time to enable application of the surface 
annual weather cycle model developed by Mitchell, (40), and (2) to establish time periods 
used in the procedures to calculate suction profiles for each case at times required to make 
calculations for the post construction algorithms. Appendix A-107 is broken down into two 
basic steps. First, the trmes required for input are entered. Second, three time periods are 
calculated using the driest month of the year as the baseline. The time period used to 
calculate the suction profile under the slab at the time of construction 1s represented by the 
variable t.onsvarys aS calculated in Appendix A-107. The time period between the time of 
construction and the time at which the soil volume change 1s sought 1s represented by the 
variable ty jconsts aS Calculated in Appendix A-107. Finally, the variable t,,,,, describes the 
time transpired from the driest month of the year to the month the soil volume change 1s 


sought. 
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All soil suction profiles with respect to surface conditions, time and depth are 
calculated using the algorithms in Appendix A-109 and Appendix A-110. The process by 
which post construction soil suction profiles are generated consist of the following basic 
steps: 

1. All values required to make calculations are retrieved for each particular moisture 
effect case and time the suction profile is sought, Appendix A-109. 

2. Suction values, for each particular moisture effect case, are calculated for 5 centimeter 
thick incremental soil layers. 

3. For every incremental soil layer cycle, variables are re-initialized and pertinent 
equations are activated based on: the location within the soil profile, the moisture effect 
case type, and the type and size of design effect implemented (If a horizontal barrier is 
implemented or if a moisture effect case 1s located a distance from the edge, Appendix 
A-104 will be used to augment the procedures described in Appendix A-110.) 

Once all soil suction profiles are generated, the soil suction profile at the slab center 
is damped from time t,,,,, to time t,, by applying the algorithms described by Appendix A- 
111. This process takes the soil suction profile under the slab at the time of construction, 
hgtab(Zistconsvidry)s and dampens it based on equations (59) and (61). These equations model 
the soil suction profile under the center of the slab at time t,,,,, so that the soil suction in 
each incremental soil layer approaches the equilibrium suction profile. The rate of 
approach depends on unsaturated soil profile properties, surface weather conditions, and 
time. The process in Appendix A-111 is broken down into two basic steps. First, all 
required variables are retrieved. Second, the dampened suction value for every incremental 
soil layer is calculated for the time transpired from the time of construction to the time the 


soil volume change is sought. 


Vertical Differential Soil Movement 
Once all soil suction profiles are calculated, the expansive volume strain theory 
developed by Lytton is applied. All moisture effect cases are applied to the same vertical 


differential soil movement algorithms which are broken down into a four step process. 
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First, all variables and quantities required for calculations are retrieved, as described in 
Appendix A-113, depending on which moisture effect case and design effect scenario is 
chosen. Second, vertical differential soil movement is calculated for soil at the edge of the 
foundation by way of steps described in Appendix A-114. Third, vertical differential soil 
movement is calculated for soil under the center of the foundation by way of steps 
described in Appendix A-115. Finally, the two quantities of vertical differential soil 
movement are summed together to arrive at the total differential soil movement for the soil 
under the slab at a specified time after construction, as described in Appendix A-115. 
Appendix B-128 is a program design layout sheet that illustrates the results of the 
post construction procedures when applied to calculate the worst case center lift mode for 
case “1-6”. The moisture effect case considered in Appendix B-128 is grass at the surface 
with a 4 ft vertical barrier, case “1-6”. The vertical soil movement calculated at the edge 
moisture variation distance under the slab is Az,,,, = -0.13 inches (shrink). The vertical soil 


movement calculated at the edge of the slab is Az,,,, = -0.61 inches. The total vertical 


edge 
differential soil movement for this moisture effect case using Table 2 as the soil profile is 
AZ,otaj = 0.48 inches (shrink). This quantity, Az,,,.,, 18 the difference between the values of 


AZ. tap ANA AZ, age: 
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MEASURED SUCTION PROFILES 


When cases exist that can not be modeled by one of the moisture effect cases 
discussed in the previous procedures, vertical differential soil movement can be calculated 
by the measured suction profile method. This is a method to solve for the differential soil 
movement by measuring the initial and final suction profiles and applying those profiles to 
the volume strain theory developed by Lytton. The procedures are broken down into three 
broad steps: (1) The measured suction values for a sample are entered and stored for later 
use, Appendix A-87. (2) The entered measured suction values for a sample are converted 
into a measured suction profile, Appendix A-101. (3) The measured suction profiles are 
used to solve for vertical differential soil movement for a specific column of soil, 
Appendix A-105 and Appendix A-106. 

The first step to applying this methodology is to input measured suction values for 
a sample of soil. There are a number of methods to measure the matrix suction in a soil 
sample. One of the more accurate means for measuring total soil suction is the 
thermocouple psychrometer. The thermocouple psychrometer sensor and data logging 
equipment is not inexpensive and requires some expertise and maintenance to operate. For 
these reasons the filter paper is becoming a more widely accepted method in practice. The 
filter paper is inexpensive, requires practically no maintenance, and requires very little 
expertise to obtain results nearly as accurate as the thermocouple psychrometer (McKeen 
1981). 

The process of inputting measured suction values occurs as part of the process to 
input moisture effect cases. Appendix B-121 is a program design layout sheet that 
illustrates the choices for the moisture effect case input. The program user would select 
‘12. Measured Suction Profile” and the next screen that appears is illustrated by Appendix 
B-122. Appendix B-122 and 123 are a program design layout sheets that illustrate 
windows environment screens for inputting the measured suction profiles. Appendix A-87 
is the process flow describing measured suction profile input. It is important to note that a 


surface suction value must be entered or a default value equal to the equilibrium suction 
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will be assumed. 

Appendix A-101 is the process flow describing the steps to generate a measured 
suction profile based on the entered suction values. The procedures are broken down into 
two basic steps: (1) The slope of a line connecting adjacent suction values within a 
sample is calculated. (2) The measured suction profile is generated, using this previously 
calculated slope, by connecting a line to adjacent suction values in the sample. Table 3 is 
an example of values used to generate a measured suction profile. Fig. 33 is a graph of the 


initial and final measured suction profiles generated by the values listed in Table 3 


Table 3 Measured Suction Profile Values 













Input Suction Values 













Initial Final 
(ft) (cm) (pF) (cm) (ft) (cm) (pF) (cm) 
Depth Depth | Suction Suction Depth Depth | Suction Suction 









3.f 9011.872 
$3.57 3715.352 


3.5 3162.278 
3.57 3715.352 







HH 6-0.00219 0.005468 


Ky -0.00109 0.001094 
iP} -0.00246 0.00164 Ko wet) = 
(tJ -0.00107 0.000574 f tany = 
20 fan = 







Applying the measured suction profiles to predict vertical differential soil 


movement is probably the simplest of all procedures contained herein. The procedures are 
broken down into two basic steps: (1) select an initial and final suction profile to 


accurately model the moisture change scenario being considered, Appendix A-105, and (2) 
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apply the algorithms to calculate vertical differential soil movement based on these 


measured suction profiles, Appendix A-106. 


Suction (pF) 
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Fig. 33. Initial and Final Measured Suction Profiles Described in Table 4 


The challenge 1s in deciding what suction profiles are used as the initial and final 
suction profile to accurately model the soil conditions which produce reasonable limit 
conditions used for design. The answer should come from knowledge about history of the 
site, to include: surface drainage conditions, soil characteristics in the soil profile, surface 
weather conditions, and type of structure placed. For example, consider the construction 
site described by the following: the site is located in College Station, the site is surface 
drained, the soil profile is described by Table 1, and the measured suction profile 1s located 
where a large oak tree was removed as a result of pre-construction clearing. The measured 
suction profile at this location 1s a dry profile. This profile will tend toward equilibrium 


and may approach the wet limit state subsequent to a wet season. A conservative approach 





2 


for calculating vertical differential soil movement is to set the initial suction profile equal 
to the measured suction profile and set the final suction profile equal to the wet limit state 
for the surface conditions that exist after the construction project is complete. This 
scenario produces heave values of vertical differential soil movement. Table 3 and Fig. 33 
are values used to obtain results for vertical differential soil movement in the example 
described above. The vertical differential soil movement in this case is 5.48 inches of 


heave. 
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DISCUSSION OF RESULTS 


The procedures to predict vertical differential soil movement were developed based 
on the breakthroughs mentioned earlier. These procedures consider moisture effect cases 
that are common with light commercial and residential structures, such as: (1) bare soil at 
the surface, (2) grass at the surface, (3) tree at the surface, and (4) flower bed at the surface. 
Additionally, these procedures include calculating the effects of differential soil movement 
caused by the introduction of design effects, such as, vertical barriers and horizontal 
barriers. The methodology to handle these moisture effect cases and design effects are 
quite delineated in detail in the process flow of detailed algorithms, Appendix A. A 
graphical representation of the program design layout is illustrated in Appendix B. 

The algorithms in Appendix A were programmed in a spreadsheet software 
package. The author used this package to analyze the moisture effect cases illustrated in 
Appendix C which consists of the soil profile described in Table 1. Table 4 is a summary 
of vertical differential soil movements for the site described by Appendix C. This 
summary illustrates the worst cases of edge lift and center lift distortion modes for all 
moisture effect cases. Post equilibrium and post construction approaches were used to 
solve for all cases so that the worst possible case of vertical differential soil movement is 
identified. Additional analyses were performed using different sizes of vertical and 
horizontal barriers to define the effects that these design features have on soil movement. 
An engineer may use the maximum cases of edge lift and center lift caused by the 
combination of moisture effect cases and the desired design effects to design an adequate 
foundation. 

Table 4 is set up in chronological order beginning with the post construction 
differential soil movements and concluding with the post equilibrium soil movements. 

The post construction numbers in Table 4 were generated for June as the driest 
month of the year and for a slab built in January (the wettest month of the year). The slab 
is placed on “wet” soil and this soil is allowed to approach an equilibrium moisture content 


based on the soil characteristics and the post construction theory described earlier. The soil 
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at the edge of the slab continues to fluctuate with the seasonal weather changes. Applying 
this post construction scenario to the “1-7 Tree” case produces the maximum differential 
soil movement compared to all other cases which cause a center lift distortion mode. The 
differential soil movement in this case is 4.32 inches and occurs the first August after 
construction. Note that a reduction to 2.44 inches of differential soil movement is made 
possible by introducing a four foot vertical barrier at the edge of the slab where the tree is 


located. 


Table 4. Summary of Vertical Differential Soil Movements for Appendix C 


Post Construction 


Driest Month of the June 
MOAN. . visissieses 
Month of January 


Construction............. 


Worst Case Center Lift Distortion Mode 


Moisture Effect Case/Design Scenario Differential Soil Movement Number of 
(units in inches) Months after 
Const. Worst 
Case Occurs 


1-4 Grass 2.34 6 
1-5 Grass with 4 ft Horizontal Barrier 1.19 6 
1-6 Grass with 4 ft Vertical Barrier 0.48 7 
1-7 Tree 4.32 7 
1-8 Tree with 4 ft Vertical Barrier 2.44 7 
1-9 Flower Bed 3.30 6 


Worst Case Edge Lift Distortion Mode 


Moisture Effect Case/Design Scenario Differential Soil Movement Number of 
(units in inches) Months after 
Const. Worst 
Case Occurs 
1-4 Grass 0.43 11 
1-5 Grass with 4 ft Horizontal Barrier 0.31 11 
1-6 Grass with 4 ft Vertical Barrier 0.23 14 
1-7 Tree 0.43 11 
1-8 Tree with 4 ft Vertical Barrier 0.23 14 
1-9 Flower Bed 0.92 12 


1-11 Flower Bed with 4 ft Vertical Barrier 0.20 12 
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Table 4. Continued 


Post Equilibrium 


Worst Case Differential Soil Movement 


Moisture Effect Case/Design Scenario Edge Heave (in) Edge Shrink 

(in) 
1-4 Grass 1.65 1.03 
1-5 Grass with 4 ft Horizontal Barrier 0.47 0.56 
1-6 Grass with 4 ft Vertical Barrier 0.07 0.09 
1-7 Tree 1.64 3.67 
1-8 Tree with 4 ft Vertical Barrier 0.07 2.36 
1-9 Flower Bed 2.89 1.86 
1-11 Flower Bed with 4 ft Vertical Barrier 0.68 0.55 
1-12 Measured Suction Profile (Point Water Source) 4.40 NA 
2-12 Measured Suction Profile (Point Water Source) 4.40 NA 
3-12 Measured Suction Profile (Tree Removal) 1.62 NA 


An inspection of Table 4 shows that the differential soil movement is greatest when 
there are no vertical or horizontal barriers in place. In the “1-4 Grass” case, the edge 
heave is 1.65 inches and edge shrink is 1.03 inches. When a four foot horizontal barrier is 
added the edge heave is reduced to 0.47 inches and the edge shrink is reduced to 0.56 
inches, as described in the Table 4 case “1-5 Grass with 4 ft Horizontal Barrier’. A further 
reduction in soil movement is obtained by adding a four foot vertical barrier as shown in 
Table 4 case “1-6 Grass with 4 ft Vertical Barrier”. This latter case reduces the edge heave 
to 0.07 inches and the edge shrink to 0.09 inches. 

As shown in Table 4, the worst case edge heave occurs in the post equilibrium case 
of the measured suction profiles produced at the locations of the leaking pipes. The 
leaking pipes produced a “wet” spot that caused the soil to heave 4.40 inches at those two 
locations. This type of heave at the edge caused the worst case edge lift distortion mode. 
Refer to Appendix B-130 for a conceptual illustration of the post construction case of 
vertical differential soil movement that can be predicted by applying the procedures 
contained herein. These procedures for analyzing vertical differential soil movements 
should be applied to a soil profile to enable an engineer to maximize resources ensuring 


that the best foundation is designed. 
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SUMMARY AND SUGGESTIONS FOR RESEARCH AND DEVELOPMENT 


The theory to predict differential soil movement of expansive soils has evolved 
since the 1970’s. The engineering profession 1s at a juncture today in which engineers 
have the technology to apply this information to maximize resources and produce high 
performing ground structures. The evolution of this theory began in the early 1970’s with 
the advent of the volume strain theory for expansive soils which was developed by Lytton 
(1973). In the decade of the 1970’s, great strides were made by Pearring (1968), Holt 
(1969), and Mojekwu (1979) to enable engineers to identify predominant clays. In 1981, 
McKeen used these newfound procedures to identify predominant clays and developed a 
correlation to link groups of clays to suction compression indices. The suction 
compression index 1s used in the soil volume strain equations developed by Lytton to 
predict vertical differential soil movement for particular types of expansive soils. In the 
1980’s, although the volume strain theory was well based, its application was limited by 
the lack of knowledge and methodology to predict the equilibrium suction and soil suction 
profiles for a particular soil profile and location. This is when work by Mitchell and later 
Gay opened the door for major advancements in this area of expansive soils. In 1980, 
Mitchell developed and applied simple mathematical methods for predicting soil suction 
profiles. Then in 1994, Gay developed a way to predict the mean volumetric water content 
for soils dependent upon the location and climatic conditions. Jayatilaka, Gay, Lytton, and 
Wray (1992) then developed and calibrated to field data a finite element method for 
determining the edge moisture variation distance. They also devised methods to predict a 
number of unsaturated soil properties using simple index properties. 

Based on the wealth of information and breakthroughs in the areas of expansive 
soils listed above, Lytton and the author made five critical advancements necessary to 
effectively and easily apply this theory: 

1. The means to calculate the equilibrium suction for a particular soil profile and location. 
2. The means to calculate depth to constant suction for a soil profile with multiple layers. 


3. The equation to calculate horizontal velocity of water flow in unsaturated soils. 





ed 


4. A method to consider the transient case for calculating vertical differential soil 
movement for expansive soils immediately after construction and before the soil under 
the center of the slab has reached an equilibrium moisture content. 

5. Procedures to apply all of the theory discussed herein to predict differential soil 
movement for expansive soils considering all of the moisture effect cases and design 
effects that have been mentioned. 

The theory discussed herein represents physical behavioral characteristics of 
expansive and unsaturated soils and 1s well documented. Those contributing to this 
vertical differential soil movement theory have approached the problem with a 
theoretically sound yet practical mindset. Additionally, this expansive soil theory has been 
developed in a way to enable practicing engineers to run simple, common and economical 
geotechnical engineering lab tests to adequately define all necessary soils profile 
information. This soil profile information is applied in the vertical differential soil 
movement theory and should be used to make sound geotechnical engineering designs. 

Further development needs to be made in the areas of enhancements to the present 
programs, materials characterization testing, in situ testing, and foundation structural 
analysis models for stiffened plates and drilled shafts in expansive soils. The program 
enhancements should include the choice to use the Gardner, as well as the Mitchell, 
unsaturated permeability options. Materials characterization testing developments that are 
needed are more convenient methods of measuring the volume change, diffusion, and soil- 
water characteristics of expansive soils directly. In situ testing developments that are 
needed include methods of measuring suction profiles on site and preferably down hole; 
lateral earth pressure coefficient; and of detecting the presence of roots and voids, either 
air-filled or water-conducting; in the field and non-destructively. Stiffened plate models 
which incorporate non-rectangular geometry need to be coupled with this program or its 
enhancements to determine design quantities of moment, shear, and deflection. Drilled 
shaft models with interface elements need to be coupled with the enhanced program to 


determine the design quantities for both shrinking and expanding soils. 
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APPENDIX A 


PROCEDURES TO CALCULATE 
VERTICAL DIFFERENTIAL SOIL MOVEMENT 
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Start of Program 
and 
Project Input 
© by Donald Naiser, Jr. 5 Aug 97 


<Ooen — 


(S$ 








Choose Program 
Function 

1. Input 

2. Make Calculations 

3. Go to Start 











3 
go to 
Pg 82 
Choose Input Functio : 
1. Project Information 7 to 
2. Foundation Layout “ Pg 83 
3. Soil Property Input 
4. Moisture Effect Input 3 5 
5. Go to Start go to 








Note to Programmer: Most of the 
input should draw from 
databases or spreadsheets in a 
“Combo Box" format to provide a 
user friendly layout. 







Project Information Input ; 
1. Input Project Name go to 
2. Input Project Date pg 82 


3. Input Project Engineer 
4. Input Site Location Address 









Note to Programmer: This 
information should be stored in a 
database file of "Projects" in a 

spreadsheet type format. This 

"Projects" file could then be used 
to organize all projects analyzed 
by this software package. 
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Foundation Layout Input General Note 


Note to user:A three dimensional rectangular (x,y,z) 

coordinate system is used to define foundation layout. The Y 
direction is defined by an angle wrt true north. The z direction is 
depth of soil (increasing positive below surface). Every part of the 
foundation structure must be in the positive (x,y) plane and the slab 
shape must be rectangular. The origin shall be defined by latitude, 
longitude, and elevation. Refer to Appendix B-118, 125, and 126 
for graphical representation. 








Specific Notes 


Inputs forlatitude and longitude 
are intended for data collection 
purposes and to link the site to site 
specific data such as; Thornthwaite 
Moisture Input data, and edge 
moisture variation data. 
Refer to the example below. 





— 


Input the slab origin 

(typically the west- | -—-—-—- - —_~ -~—-. 
most and the south- 

most corner of the 

Slab. 








—— <= The origin elevation is simply used as 
Input the elevation of a reference for 3D computations. 
the origin. Refer to the example below. 













This hasno bearing on calculations. 
This simply allows the user to establish 
building location/alignment. The value 
should be entered clockwise from true 
north. Refer to the example below. 


Input the "y" axis 
direction 


Input the Local TMI 
TMI = 


Hie 


Input the number of — _— . 
Local weather cycles/y Example: Cut-out 1 from Appendix B-118 


n= 


2 ie Specific Note 

Input the information me | 

for the corners of the ae Alt cells in the "corner information" 

Slab. spreadsheet are critical. The user 
must not be allowed to continue without 
first entering this information. 





go to 
Pg 82 
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Soil Layer Property 
Input 


(Boring sample 
Input: Number of boring samples counter) 
B= b= 


b=b+1 
s = QO (layer counter) 


input location of boring, 
number of soil layers, and 


abut Soi! Layer Properties: depth of boring sample. 
____ Depth to next layer(ft) 
, =__Liquid Limit forlayer s (%) 
= Plastic Limit forlayer s (%) vee 
___ Percent passing 2 micron (%) ___ (number of layers) 
|___ Percent passing #200 sieve (%) BD, =~ (depth of boring 
Yds=— :_ Dry density (pcf) sample) 
Refer to Appendix B-120 for program Refer to Appendix B-119 for 
design layout. program design layout. 

















Calculate additional 

soil property quantities 
required to make 
additional 

computations. 








Options: 

1. Go to Soil 
properties 
display 

2. Go to Start 















Give user the option to view 
spreadsheet of other calculated 
material properties used in volume 
strain, suction and water flow 
velocity calculations. Refer to 
Appendix B-126. 











7 
[ 
eae 

anes 


1. Bare soilat surface 

2. Bare silwithhorizontalbarmer = 
S. Bare soil with vertical barrier 
4. Grassat suriace 

5. “Gris with honizontalbarmer 


6. Gn ss with vertical barnier 


7 Tree near foundation edge 


: : } : 
drmrvem net and nen nd rd trannies mpadsmnrensrosredrerns rpwennn haces nrAeeecessndinnt see By 


Dansnndie 


g. Tree wath vertical barrier 
|S. Flower bed near foundation edge — se 
10. Flower bed with horizontal barrier 

11. Flower bed with vertical barner 
12. me mered 5 suction apes 5 


Se riple: Bitcit 2 from Appondls B-121 


Program Design Layout 





GRAPHIC DISPLAY 


——— 
Conan 


dence EL 


i 
i 
3 
i 
i 
i 
i 


een es i eee aD x cece ees oneonnt st care essere 


Example: Cut-out 4 from Appendix B-1 21 









: 
See erection. 
t 


Nertical 
Barriec Ps 
(Runs along the full | 
tength pf the Rastedge of F 
the foundation.) Bi 


Moisture Effect Input 















Select Moisture 
Case: Refer to Cut-out 

2 below for possible 
selections, or 
return to Start. 









As the Moisture Effect Case type 
are selected a database 
spreadsheet row is generated an 
a counter is initiated for each cas 











Example 













12 start 
40 Start For example, the first "flowerbed" moisture 
go to go to effect case selected would be given the 
pg 87 pg 82 “case number’ name (1-6) [and the second 






would be (2-6) and so on]. A row would be 
established to contain all pertinent cells of 
information required for computations made 
by this software package for that particular 
case. Refer to cut-out 7 on the next page. 












Prompt user to input a "start" and "en 
location for all moisture effect cases. 





von 


Example 






Example 





For example, the "start" and "end" points for 
the flower bed are (40,0) and (0,50) 
respectively. If the moisture effect is a point 
source then the two inputs will be the same. 
Refer to the cut-out 3 and cut-out 4 for 
location of flowerbed on grid and 

example of coordinate inputs respectively. 








Example 






go to 
pg 86 


85 





86 


Moisture Effect Input 
Continued 


papoanspsoosseespensecasenaanastianienses 


nm PVUrl 





1, Bare soil st rurtace 
«Bare sil with dorizontal barrier, 





2 
=] 5. Pare soil with vertical barrier 
4. Grassat suriace 
75. Grass wath horizoatal barrier 
6 
? 
9 






4416. Gross with vertical barrier 
] 7. Tree near foundation edge 


Perform the following 







function for moisture effect 
cases 7 thru 11. 


Tree wath vertical barner in 





gi | 9. Flower bed near foundation edge 

| 20. Flower bed with horizontal barrier 
11. Flower bed with vertical barrier 
12. Dea sured suction profiles 















Example: Cut-out 5 from Appendix B-121 










—— 
—— 


Prompt user to input the depth of the 
a moisture effect zone. 
Example d 





For example, the depth of the flower bed in 
this case is 4 feet. Refer to cut-out 5 shown 
above. 







zone ~ 








Refer to cut-out 6 for example of edge ——_———-Example_ 
location input for the flower bed. ; 






Is the effect 
case at the edge 
of slab? 









——~ 


Prompt user to enter 
location wrt slab. 
(-) number indicates 
the zone is under the 
slab, (+) number 
indicates the zone is 
away from the slab. 
d 













When yes is selected 
this indicates the effec 
is at the edge of the 
slab. 











source _ 








Example: Cut-out 6 from Appendix B-121 








Provide the user opportunity to 
review spreadsheet of moisture 
effect input data. Refer to cut- 






Does 
the location have 
a barrier? 



















es Example/ 
Note 
No 
Default size of barrier 
is zero. Input size of barrier. 
dwar = O (vertical) i 
dibar = O (horizontal) divar = 
6 
go to 
g 85 
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Measured Suction 
Profile Input 


10 
from 


pg 85 


Input number of Suction 
Measurements (including 
the surface value) 

N= 


Initialize Suction 
Measurement Counter: 
c=1 












Note: Measurement one (h,) 
must always be at the surface and 
D, is always 0. If no value is 
entered a default value ofh,, 
(equilibrium suction, calculated 
on pg 90) will be assumed. 


Input Suction Value 














Input the Depth to 
Measurement "C". 
D. = (cm) 


Input the suction value 
@ Measurement "C”. 
he = (pF) 
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Soil Property a 





eer 
Calculations Initialize Variables | 
M, (retrieved from pg 84) 
i = 0 (increment counter) 
increment = 5 (cm) 
is s = 0 (layer counter) 
po e4 z,= 0 (depth to the i increment) 
=r > a « €¢(Z) =0 
Ene(Z) = 0 
Input Quantites 








(Retieve values from pg 3) 



















Bey LS fied 

PL(z) = Ee z= z,, + increment _ 
H(Z,) = 1, oe 
u(iZyia, 










Ya(Z) = Yas 
eee 


Activity Ratio Cation Exchange Capacity (in percent) 


_PI(z,) 
Clay(z; ) 









Plasticity ines an percent) 
Plz) = LL(z) - PL(z) 


Percent Fine Clay (in percent) 
CI u(Z,) 
ay(z;) = ~~ 
VA Z. 
Slope of the Suction vs Gravimetric 
Water Content wet Line 








The Larger of the CEC(z )= = PL(z,)'. 17 
two quantities CEC(z j= LL(z, jo 912 


Ac(z;) = 













Cation Exchange Activity 


CEC(z, 
CEAc(z,) = aa 










Matric Suction Compression Index 








Use calculated values of Ac(z , ) and 
CEAc(z,) and read values of y,,., from the 
chart in Figure 15 developed by McKeen. 


















¥n(Z;,) = Clay(z;)-¥ soo 





Yes 
Mitchell's Diffusion Coefficient (cm 2/s) 
a(z,) = 0.0029 —0,000162-S(z,)-0.0122-y,(z,)) 


Edge Moisture Variation Distance 
Center Lift 


Mitchell's Unsaturated Permeability 
a(Z;)-¥ 4(Z,) 












h,| = 0.4343- 





Yes 




















Edge Moisture Variation Distance 
Edge Lift 

















Use TMI (input from pg 83) and a{z)) 
(calculated above) and go to Figure 16 
to find the appropriate edge moisture 
variation distance for center lift [e ,,,(Z)]. 


Use TMI (input from pg 83) and a(z,) 
(calculated above) and go to Figure 17 
to find the appropriate edge moisture 
variation distance for edge lift [e _.(Z)]. 


“> 
IS Cn_(2) > ee 
Yes 
“| One BIne(Z) 





















Cnc (2) =Enc (2, ) One (Z)=En (Z,, ) 









IS Cnc (2) > @ng(2,4) yey 


—_ Yes 
ial Png (2) = Ene (2) 








2 (zeae 
il Oo 
a vs 
go to No—<= is s <=M, a oe 





pg 84 
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Calculations for Calculations for 
Post Equilibrium Post Construction 
Moisture Effect Moisture Effect 
Cases Cases 


















1. Post Equilibrium; 
2. Post Construction, 
3. Go to start. 














Calculate Equilibrium 
Suction for Bare Soil 


Select 3 


Input Construction 
Time Constraints and 









and Soil with 
Vegetation @ the a Time 
surface. aselines 










Calculate Depth to 
Constant Suction(Zm) 
for all Moisture Effect 
Cases. 










Calculate the Suction 
Profiles for all Moisture 
Effect Cases. 









Calculate the 
Differential Soil 
Movement for all 

Moisture Effect Cases. 






Calculate Suction 
Profiles. 












Calculate the 
Differential Soil 
Movement for all 

Moisture Effect Cases. 








13 
go to 
Dg 89 





Equilibrium Suction 
h 


m 


Input: 
1. Inputd,,, (units of cm), 
2. Inputh,, (units of pF), 
3. Input ary (units of pF). 


If input values are not known, then 
the following default values will be 
used; 

dm = 30 cm 

h,, = 2.5pF 
hy, = 6.0pF (bare soil) 
Navy = 4.5pF (vegetation @ surface) 

6, = 0.55 

6,= 0.04 

A = 160 

B = 0.68 
if the USC is known, then go to Table 
2 to find the corresponding values of 
6,, 6, A, and B. 


h, =h,+.01 
h,, =h,, - .01 ws 
No 


No 





Yes 


Calculateé,, & @',: 





90 


Retrieve TMI from 
pg 83 


Calculate 6,,: 


6 teal g 
A 


Fr 








Calculate Bary: 







abel 
se) Aes 


B 
pe bol 


Calculate constants: 
T = TMI + 60 
y = 0.039337d_, + 1.357033 
d, = 0.449079d,, + 0.304560 
T, = 0.062651d, + 59.53593 





Initialize h, & 8: 
he = (ngs hy y/2 
5=1 








Depth to Constant 
Suction 


zm 







The algorithms used to compute the 
depth to constant suction may vary 
depending on the moisture effect 

case being considered. 


Go to the proper depth 
to constant suction 
algorithm depending 
on which moisture 
effect cases are 
entered. 

















Calculate z,, 
for the following cases 
(if the case exists): 

1. Bare soil 

2. Grass at surface 

3. Tree 

4. Flower bed 

5. Measured suction values 
6. All Z_s calculated. 


20 
go to 
pg 89 


3 


9] 











Depth to Moisture Active 


Depth to Constant 
Suction 
For Bare Soil 


15 
from 


pg 91 


Define Variables 
h,,, (bare soil) -h,, (pF) 
layer number considered given an 
initial value of zero. s=O 
the vertical thickness oflayer s (cm) 
depth to constant suction given an 
initial value of Zero. 2 =0 
depth to layer s (retrieved from 


pg 84) 








Solve for Z, 


nal 


2 
Z final =—In ue : 


o 






Zone 


2m = Zm i Zrinal 











oe 


Depth to Constant 
Suction 
For Grass @ 
Surface 


16 
from 


pg 91 


Define Variables 
h,, (vegetation) -h,, (pF) 
layer number considered given an 
initial value of zero. s=0 
the vertical thichness oflayer s (cm) 
depth to constant suction given an 
initial value of Zero. 2 =0 
depth to layer s (retrieved from 


pg 84) 


Solve for P.(the ratio of the suction 
difference between the top 
increment and bottom increment fo 
layer s). 


n = Weather cycles per second. 
a, = Mitchell's diffusivity for layer s. 


Calculate U,. The 
difference between the 
maximum and 


minimum suction 


values at depth L. 
Ulu; 








YS 


Depth to Constant Depth to Constant 


Suction Suction 
For Flower Bed @ 
rs For Tree @ Surface or Surf 
from from urTace 
pg 91 pg 91 


Define Variables 

Naty (PF) 
layer number considered given an 
initial value of zero. s=0 
the vertical dimension of layer s(cm) 
depth to constant suction given an 
initial value of zero. z_,=0 

= input value for depth of flower bed 
zone (retrieved from pg 86) 
depth to layer s (retrieved from 
pg 84) 


Define Variables 
hy (vegetation) =h., Wtpr) 
layer number considered given an 
initial value of zero. s=0 
the vertical dimension of layer s (cm) 


depth to constant suction given an 
initial value of Zero. z_.=0 

input value for depth of root zone 
(retrieved from pg 86) 

depth to layer s (retrieved from 
pg 84) 





Yes 










Is on = done = ZmtZs 


= ee N 
Zz, Zmt2s Sad = 


Solve for P, (the ratio of the suction 
difference between the top and 
bottom increments for layer s). 







zZ_.=d 


m zone 







n = Weather cycles per second. 


Solve for Z,,4 a, = Mitchell's diffusivity for layer s. 


Al I 
‘ae al a 


U n-1 











—Z, 


Pze V% 














Calculate U,,. The 
difference between the 
maximum and 
minimum suction 
values at depth L,. 
Wir a if ; GF 


Depth to Moisture Active 
Zone 
Zn = 2m i Z final 





94 


Depth to Constant 
Suction 
For a Measured 
Profile 









Define Variables 







layer number considered given an initial value of zero. s=0 
depth to layer s (retrieved from pg 84) 





” 










Ze depth ofincrement s (in cm), 

Z,= depth to constant suction given an initial value of zero. z,=0, 

h= the c' suction value for a sample (retrieved from pg 87) 

N= Number of suction measurements per sample (retrieved from pg 87) 
D,= depth to the N ™ suction measurement (retrieved from pg 87), 
h,,(D,)=equilibrium suction at depth I) (retrieved from pg 96) 





( 
= _ absolute value of (h,, - h_,(D,)) 


t= 


= 
= 







Z,=2m Value for the surface 
conditions most similar to 
the surface conditions for 
the measured profile but 
not less than z, for bare 

soil. 












is Z,>=Z,, value 
of the surface conditions mos 
similar to the surface conditions 
for the measured profile but no 
less than z, for bare soil? 






No 








Z= Lo 2. 


Solve for P, (the ratio of the suction 
difference between the top and bottom of 
Z=L-L, soil increment s) 


n = Weather cycles per second. 
a, = Mitchell's diffusivity for increment s. 


nn 
~—Zz oo 
"Va, 


P.=e 


Yes Calculate U,. The 


difference between the 
Solve for Z,_. 
0.2 l 
Z final — o(*2] . nt 











maximum and 
minimum suction 
values at depth kL. 


Uae 





Depth to Moisture Active 
Zone 


Zn a Zm a 2Zrinal 











oS 


Calculate Suction 
Profiles 
Menu Page (2) 


The algorithm used to calculate a 
suction profile wrt depth varies 
depending on the moisture effect 
case and surface conditions. 







Go to the Propersuction profile 
algorithm depending on which moisture 
effect case is being considered. 





















Calculate 

Suction Profiles 

for the following cases 
(if the case exists): 

1. Equilibrium (for all cases), 

2. Dry Limit Suction Profile 
(for bare soil), 





3. Dry Limit Suction Profile 3 go to 
(for grass @ surface), pg 98 
4. Dry Limit Suction Profile 4 
(for tree @ surface), 
5. Dry Limit Suction Profile 26 
(for flower bed @ surface), 9 eee 
6. Wet Limit Suction Profile 2 
(for bare soil), 6 


7. Wet Limit Suction Profile 
(for grass @ surface), 

8. Wet Limit Suction Profile 
(for tree @ surface), 

9. Wet Limit Suction Profile 
(for flower bed @ surface) 

10.Measured Suction Profiles 

11.All Suction Profiles 

calculated. 


32 
go to 11 
pg 89 








Equilibrium Suction 
Profiles 











Surface 
condition: 
1. Bare soil, 
2. Grass, 

3. Tree, 

4. Flower Bed 


Tree 
@ Surface 
Retrieve: 

h,, and z,, previously 
calculated by the 
processes shown on pg 
90 & 93 respectively. 


Grass 
@ Surface 
Retrieve: 
h,, and z,, previously 
calculated by the 
processes shown on pg 
90 & 92 respectively. 


Bare Soil 
@ Surface 
Retrieve: 

h,, and z,, previously 
calculated by the 
processes shown on pg 
90 & 92 respectively. 














Convert: 
Suction in pF to 
centimeters. 


The process to calculate 
h,,(Z,) is the same for each 
case. The only difference 
is the h, @ the surface 


and the z, for each case. h_ (cm) = 10%*F) 


No 


21 
Goto Yes Z=Z, 
20 95 Ny (2) = Pg 





96 









Flower Bed 
@ Surface 
Retrieve: 
h,, and z,, previously 
calculated by the 
processes shown on pg 
90 & 93 respectively. 


Initialize Variables: 


i=0 
Z, = 0 (cm) 


increment = 5 (cm) 


Nn Zo) = ne 


ian 
, + increment 
(Z,,) - increment 








Dry Limit Suction 
Profile 
For Bare Soil & 
Grass @ Surface 


Bare Soil 
@ Surface. 
Retrieve Values Calculated in 


previous process flow sheets: 


Z, from pg 92, 
h,, from pg 90, 
h,,(Z,) from pg 96, 
n from pg 83, 
a(Z,) from pg 88. 





oe 


Grass 
@ Surface. 
Retrieve Values Calculated in 
previous process flow sheets: 


Z, from pg 92, 
h,, from pg 90, 
h,(Z,) from pg 96, 
n from pg 83, 
a(Z,) from pg 88. 


Input Dry Suction Value @ the surface. 


an (pF) 


max = Ngy = 6.OpF (bare soil default) 
max = Nyy = 4.5pF (vegetation @ surface default) 


Initialize Variables: 


i = 0 (Increment counter) 
increment = 5 (cm) 
z,= 0 (Depth to I" increment) 
h waxtZ) = h,,, (input/default) 
s = 1 (Layer counter) 


Ziop = 9 


i=i+1 
Z, = Z,, + increment 





-( 
h,»(Z;) = h,,(Z;) + Ue 





=< 


Yes 





Calculate Maximum Suction at depth z;: 


Zi ~Ztop ) 


Retrieve 
The following Values 


M,, from pg 84 
L. from pg 84 









U, = inax(Z4) ~ Hin(Z.4) 
Zp = 2. 
s=s+1 


go to 
pg 95 
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Dry Limit Suction 
Profile 
For Tree and Flower 
Tree @ Surface. | Bed @ Surface 


Retrieve Values Calculated in 
previous process flow sheets: 


25 z,, from pg 93, - 
h., from pg 90, 7 


h_(Z) from pg 96, Retrieve 

















n from pg 83, Input Dry Suction The following Values 
a(Z,) from pg 88. Value @ the surface. 
"Pow = (pF) M, from pg 84 
: = hy = 4.5pF ls, from pg 84 
Flower Bed @ Surface. (Default). done from pg 86 










Retrieve Values Calculated in 
previous process flow sheets: 


26 z,, from pg 93, 
h,, from pg 90, 

h_(z,) from pg 96, 
n from pg 83, 

a(Z,) from pg 88. 


Initialize Variables: 


i = 0 (Increment counter) 
increment = 5 (cm) 
z,= 0 (Depth to #” increment) 





Hnax(Z) = Ana, (Input/default) 
U 


o max” ''m 


s = 1 (Layer counter) 












a se 
ea ep oe il Be 
ae p< isz>=L.? Yes 
Z, = Z,, + increment ee a 


Hnax(2)) =h 


= Pnax(Zi1) - An (Z.4) 
Ztop i Zi 
s=sti1 





es —— Yes——_< as 


No 






Calculate Maximum Suction at depth z: 


n-7 


—(2;-2 top ) Ae 
hax (Z,) = h,,(Z,) + U,€ 


a2 : ~? “21 
—___—__—_ —_—Yes- isZ,<=Z, ? No go to 
pg 95 





Wet Limit Suction 
Profile 
For Bare Soil, 
Grass, & Tree @ 
Surface 


Initialize Variables: 


i = 0 (Increment counter) _ 
increment = 5 (cm) 
z,= 0 (Depth to i ® increment) 
henin(Z) = Ayi,(Input/default) 
U, = Nm(Zo) 2 Dinin 
s = 1 (Layer counter) 


Z,=Z_, + increment 





Yes 


Bare Soil 
@ Surface. 
Retrieve Values 
Previously 
Calculated: 


z, from pg 92, 
h_, from pg 90, 
h_(z,) from pg 96, 
n from pg 83, 
a(Z,) from pg 88. 





Calculate Maximum Suction at depth z ;: 


-(z 
h_4,(Z;) = h, (a4 Use 


Grass 
@ Surface. 
Retrieve Values 


Previously 
Calculated: 


Z, from pg 92, 
h,, from pg 90, 
h_,(Z,) from pg 96, 
n from pg 83, 
a(Z,) from pg 88. 





Input Wet Suction 


Value at the surface. 


ee (0) a) 
Dinin = Ny = 2-5PF 
(default) 


Retrieve 


The following Values 


M, from pg 84 
L, from pg 84 












a 
mia en 


21 
No go to 
pg 95 





U, 7 He (Z.4) 7 Penin(Z,-1) 


oo 


Tree @ Surface. 
Retrieve Values 
Previously 
Calculated: 


z,, from pg 93, 
h,, from pg 90, 
h,(Z,) from pg 96, 
n from pg 83, 
a(Z,) from pg 88. 





Ztop aa 


s=s+1 
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Wet Limit Suction 
Profile 
For Flower Bed @ 


(0) Surface 


Flower Bed @ Surface. 
Retrieve Values Calculated in 


i Retriev 
previous process flow sheets: etrieve 


Input Wet Suction The following Values 
z,, from pg 93, Value at the surface. 
h_ from pg 90, fa (pF) M, from pg 84 


L. from pg 84 
h,,(z;) from pg 96, d,_ from pg 86 


n from pg 83, zone 
a(z;,) from pg 88. 


Initialize Variables: 


i = 0 (Increment counter) 
increment = 5 (cm) 
z,= 0 (Depth to i'" increment) 
hain (Z;) = A nin (Input/default) 
oF hn (Z,) Si Dinin 
s = 1 (Layer counter) 
Zio) = O 


Mi 










U, = hy (2.4) > Minin (2-4) 
Ztop — ie 


Ss=s+1 





i=i+1 
; Yes 
Z,=Z,_, + increment 


eZ) a lee es 


min 





Calculate Maximum Suction at depth z;: 


ner 


—(2)~-Ztop ) a(z;) 


h in(Z;) = h,,(z,) + U,e 





Yes 
pg 95 





Suction Profile 
For a Measured 
Profile 


Retrieve The Following Values 


D, = depth to measurement c (pg 87), 
Netmun = Measured suction values from c thru to 
the N,, measurement (pg 87), 

c = arbitrary measurement counter. 


Z,, = depth to constant suction (pg 94), 
(1) N = number of suction measurements (pg 87), 











Initialize Variables 


i = 0 (increment counter) 
increment = 5 (cm) 
z, = 0 (depth to ith increment) 
c = 1 (suction measurement counter) 






jsit+t 
2,=2Z,, + increment 


h(z,) = h(z,,) + 2,“Slope, 







10] 


Initialize Counter 


c=0 


Calculate Slope 
h +) h, 


c 


a D.,, ~~ D, 






Slope 





Yes 


21 
go to 
pg 95 
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Calculate Soil 
Volume change (2% ) 
Menu Page 


The algorithm used to calculate soil 

volume change varies depending 
on the moisture effect case and 

surface conditions. 






Go to the Proper soil volume change 
algorithm depending on which moisture 
effect case is considered. 















34 
go to 
Dg 103, 

Calculate 
Volume Change 
for the following 
cases (if the case exists): 1-3. 4.:6° 

1. Bare soil at surface, 7,8,9,1 

2. Bare soil with horizontal barrier, 

3. Bare soil with vertical barrier, 

4. Grass at surface, 

5. Grass with horizontal barrier, 

6. Grass with vertical barrier, 2,5, 

7. Tree at surface, 10, 12 mat 

8. Tree with vertical barrier, 

9. Flower bed at surface, 

10. Flower bed with horizontal barrier, 

11. Flower bed with vertical barrier, 

12. Moisture effect some distance from 13 

the barrier/slab edge, 
13. Measured Suction Profiles, 
14. All quantities calculated. 
36 
go to 
Dg 105 
2 
14 go to 


pg 89 
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Soil Volume Change 
for Bare Soil, Grass, Tree 
and Flower bed at the 













surface with | 
Vertical Barrier | wy < 
Solve Maximum 
. value for: 
34 1. Shrink, 
from 2. Heave. 
g 10 Z 
Retrieve Variable Initialize Variables Initialize Variables 
k_ = .33 = .67 
dear fFOM pg 86 f= 5 t's 
h(z,) = h,(z) (from pg 96) h,(z,) = h(Z,) (from pg 96) 
h(z) =h_, {z) h.(z ) = h,(Z) 
——e [ee (z) for bare ‘soil/grass (Hon (Z,) for bare soll, grass, and 
Initialize Variables go to pg 97] tree go to pg 99] 
6, = 40 (cm) [h,.... (Z,) for tree/flower bed [h__,, (Z,) for flower bed go to pg 100] 
Az= 0O(cm) go to pg 98} 


i = 0 (soil increment counter) 
z,= 0 (depth to i” increment) 


increment = 5 (cm) 
i=i+1 
Zz; = Z,, + increment 


(4¥) . -1,(2)lopy( 2 ze a 


GS 


















Convert all Suction in pF to centimeters. 


h(cm) = 10°F) 





he(2) = hy(2) 










a 
Yes 


+). =~ ¥n(Z; ogy Re 
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Soil Volume Change 
for Bare Soil, Grass, Tree and 
Flower Bed at Surface with 
Horizontal Barrier 













Retrieve Values 
Note: hoizontal barrier can no 
be used for a tree. 


35 
a dinar from pg 86 
g 10 q source from pg 86 
Default values are: 
dibar> 0 


source ~ 


Solve Maximum 
value for: 
4. Shrink, 
2. Heave. 


Initialize Variables Initialize Variables 
| = O (soil increment counter) i = 0 (soil increment counter) 
z,= 0 (depth to I increment) z,= 0 (depth to I" increment) 
increment = 5 (cm) increment = 5 (cm) 
Cnc = center lift (from pg 88) One = edge lift (from pg 88) 
| Ne i(Z) = Anax(Z) (PG 97 for grass & bare soil he a(Z) = Ny ia(Z) (PQ 99 for grass & bare soil) 
h..q(Z) = Hya,(Z,) (Pg 98 for flower bed) he a(Z) = Npin(Z,) (Pg 100 for flower bed) 
h,,(z,) = equilibrium suction (from pg 96) h,,(Z) = equilibrium suction (from pg 96) 
edge = Cnc ~ Anbar deage = Ene ~ Gnbar 
















Convert all Se 
Suction in pF to centimeters. a 4 i=ir+i1 


ee Z,=Z,, + increment 
h(cm) = 10°? ~ 













Calculate Horizontal Velocity 


a K,(Z;)- h, | 
ula) d pees 


Cnecle source 


1S Z5—7.7 No goto 
an 20? g 102 


: 


Use h,,,,(Z,) in the place ofh,(z)) 

on pg 103 to make volume change 
calculations for bare soil at the surface 
with a horizontal barrier. 
















Calculate Suction at Foundation edge 








—Vx (2; )-d edge 


Neace (Z;) = ti(Z) 6 K,(z,)-|ho| 








Soil Volume Change 
for Measured Suction Profiles 





values of 


h(z,) and h,(z,). 





No 
(Shrink) 


Initialize Variables 


2 
5 
from pg 96) 
h-(2) = Dinax(2) 
(z,) for bare soil/grass 
go to pg 97] 
a» (Z,) for tree/flower bed 
go to pg 98] 


{h 


max 


(h 


m 


1. (h, & h,)= entered 








Options: 


values, 


2. h, = equilibrium profile 


h.= entered value, 


3. h,= entered value 


h.= equilibrium profile. 


(Go to pg 87 input 





Go to pg 101 to obtain 
values of h.(Z)). 

Go to pg 96 to obtain 

values of h(z)=h, (Z). 





Initialize Variables 
G, = 40 (cm) 
Az = 0 (cm) 
i = 0 (soil increment counter) 
z,= 0 (depth to i increment) 
increment = 5 (cm) 


36 (a) 
go to 





105 







values of h(z). 
Go to pg 96 to obtain 
values of h.(Z,)=h,,(Z). 





Yes 
(Heave) 


Initialize Variables 


8 
hz) = h,,(z) (from pg 96) 
h-(Z) = Dinin(2Z,) 
(Z,) for bare soil, grass, and 
tree go to pg 99] 
(Z,) for flower bed go to pg 100) 


fh 


min 


[h 


min 








Soil Volume Change 
for Measured Suction Profiles 
(Continued) 


Yes 


33 
goto No 
Dg 10 


; AH 
Az = Az+ increment: (2H) 










isi+t 
Z,=Z,, + increment 
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. . 
. 


Convert all Suction in pF to centimeters 


h(cm) = 10°C) 


AV 5,(Z\) 


ol =-14(2,)}og4[ = 
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Post Construction Time 
Constraint Input and 
Calculations 


Input the driest month of the year a 
Input the month of construction ipaS 
Input the month soil volume change is sought t,,= 


Definition of Variables 
Driest month of the year. 
toonst= -Ehe month the slab was constructed. 
t= The month the soil volume change is sought. 





t =Time elapsed between the driest month and 
and the month of construction. 

= Time elapsed between the months of 

construction and the month the soil volume 

change is sought. 

\eary= Time elapsed between the driest month and the 

month the soil volume change is sought. 

testi1= Dummy variable. 

test2= Dummy variable. 


const/dry 






Set text values equal to 
numberical values fort,,, 


Conte and t.. 


"Jan"=0 "Jul"=6 
"Feb"=1 “Aug"=7 
"Mar"=2 "Sep"=8 
"Apr'=3 "Oct"=9 
"May"=4 “Nov"=10 
"Jun"=5 “"Dec"=11 







Nicicoret 











test2 = 12 
Ye 


t 


Ss 
const/dry = test1 + test2 N —<stestt <0? > test1 = Leonst . lary 
[ware 









test2 = 12 


test2 = 0 


Yes 





Cecawm — esti + test2 test? =f. - tone 





No 











Lay = Wercenate” “coftiuaiy 






Convert Time in 
Months to Time in 
Seconds 









Calculate Post 
Construction 
Suction Profiles 


Menu Page (8 


The algorithm used to calculate a 
suction profile wrt depth and time 
varies depending on the moisture 
effect case and surface conditions. 







Go to the Proper suction profile 


algorithm depending on which moisture 
effect case is being considered. 















Calculate 
Suction Profiles 
Based on the time entered 
on pg 107 for the following 
cases (if the case exists): 






1. Bare soil at the surface, 
2. Grass at the surface, 

3. Tree at the surface, 

4. Flower bed at the surface, 
5. Under slab profile at t.,...,, 
6. Under slab profile at t,,, 
7. All Suction Profiles 
calculated. 


47 
go to 7 
pg 89 


108 





1 pg 109 
= 42 
go to 
3 pg 109 
4 43 
go to 
5 Dg 109 








Retrieve and 
Intaialize Variables 
for Suction Profile 

Algorithms 


Bare Soil at Surface. 
Retrieve Values from previous 
process flow sheets: 


n from pg 83, 
treary (SeC) from pg 107, 
z., from pg 92, 

h. from pg 90, 


h,, from pg 97, 


h,, from pg 99, 

M,, from pg 84, 

L. from pg 84, 
tat any: 


h_ (z) from pg 96, 


n from pg 83, 
a(Z,) from pg 88, 
h, ry from pg 97, 
tz) from pg 97 
h,,, from pg 99, 
M, from pg 84, 
L, from pg 84, 

t= tary: 


Under Slab Profile at time of 
construction. 
Retrieve Values from previous 
process flow sheets: 





N onktfdry nae from pg 107, 





™~ “constidry 





Surface conditions 
under slab prior to 









Use the time variable listed 
above (t=t 








construction: 
constidry) In the place of . 
time variable (t=t,.,,,,) and 1. Bare soil, 
perform suction profi ile 2. Grass, 
calculations depending on the 


moisture effect conditions that 
existed at the surface just prior to 
placing the slab. 


3. Tree, 
4. Flower bed 


Grass at Surface. 
Retrieve Values from previous 
process flow sheets: 


treiary (SEC) from pg 107, 
zZ,, from pg 92, 
h,, from pg 90, 
h_,(Z) from pg 96, 


a(Z,) from pg 88, 


Hinax (2) from pg 97, 





45 
go to 


pg 110 
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Tree at Surface. 
Retrieve Values from previous 
process flow sheets: 


Cciary (sec) from pg 107, 
Z,, from pg 93, 
h,, from pg 90, 
h,,(Z,) from pg 96, 
n from pg 83, 
a from pg 88, 
from pg 98, 
neta from pg 98, 
h,, from pg 99, 
h.(z) from 99, 
M,, from pg 84, 


i from pg 84, 
te ty 





Flower Bed at Surface. 
Retrieve Values from previous 
process flow sheets: 


treiary (SEC) from pg 107, 
z, from pg 93, 


h. from pg 90, 
h_ (z,) from pg 96, 
n from pg 83, 
a(Z,) from pg 88, 


a’ from pg 98, 

hax (2) from pg 98, 

h,, from pg 100, 

M, from pg 84, 

L. from pg 84, 
Se (olan: 
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Suction Profile 
Algorithms 


Initialize Variables: 






i = 0 (Increment counter) 
increment = 5 (cm) f gs 
z, = 0 (Depth to ith increment) 
U,=ha,y a Dy tae 
s = 1 (Layer counter) 


. i=i+1 
Z, = Z,, + increment 


eee —— U,, =h,, (Z,, ee nin (Z,; ) 


= IE 2s DeVoe oo = Minar (1) ~ An (2.4) “top 


a 


ie cee a 





eg ott 


Sol 








For tree, if h(z,,t)<h, (z,) | { I Zaps oe en 2a | 
then U,=U,, ESD U 


0 


else all others 
No U,=U,. =t 


new 

















i 2 > Noi h(z,t)>=h,,(2i)? > 
a ee ae _ ee 


Yes “is h(z, t)>=h,, (2)? > 


1. Bare — 


2. Tree 
3. Flower Bed 












sO 
No 


J 


U, = ha, (Z,,) ¥ Pnin (Z,,) 





Yes 


Calculate Maximum Suction at zand t: 


nh: 7 
h(z;,t) =h,,(z;) + U, -co 2M Zi Ze 2) 


Calculate Maximum Suction at zand t: 









nx 





-(z, — Zion ) 


h(z,,t)=h,(z;)+ U,e a co Mt (272, 


h,,a(Z,,t)=h(z,,t) for bare soil, 
Norass(20)=h(Z,t) for grass, 

HN, ee(2Z:t)=h(z,t) for tree, 

hp, (2,t)=h(z,t) for flower bed, 
Njap(2t)=h(Z,,t) for under slab, 
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"Damped" Suction 
Profile Underneath 
the Slab 


"h Z § 
slab| ver dry) Retrieve Values from previous 
process flow sheets: 


h_,(z,) from pg 96, 
fee tee ) from pg 109/110 
a(z,) ft rom pg 88, 
n from pg 83, 
h, from pg 88 
L. le pg 84, 
teonstiary “(sec) from Pg 107, 
Ccicones (sec) from Pg 107, 
cldry (SEC) from pg 107, 


z,, from pg 92 or 93, 
i=i+1 B =1 (an exponent in from 
Zz, =Z,, + increment mitchell’s equation), 
t,(z,) time when the value U(z,) 


is reached, 
U,(Z,) is a suction value half way 
between 
h,,(z,) and h,,,.(Z 


Initialize Variables: 


i = 0 (Increment counter) 
increment = 5 (cm) 
z, = 0 (Depth to f* increment) 












i? tongtiary) 







U,(z,) = eee) = h,,(2,)] 


2 


a(Z,) 


SI 


-In(O 5 








U,(z,)-h,,(z,) ae 
a= east — a step (2 jo Cconaera) e eee | 


Mey 2 qteree eno 7 U,(z;) 0.8686 


)-a-t a 


vc/const 


ny (Z;) + A tap (Z; ’ U once 
Iga a 


vc/const 


39 
Neo go to 


Dap (Z; ? tcrany) - 





a) oe 
a —_—_ 
—_ =_ 





Post Construction 
Soil Volume Change (4 ) 
Menu Page 







The algorithm used to calculate soil 

volume change varies depending 
on the moisture effect case and 

surface conditions. 


Go to the Propersoil volume change 


algorithm depending on which moisture 
effect case is considered. 













Calculate 
Volume Change 
for the following 
cases (if the case exists): 
. Bare soil at surface, 
. Bare soil with horizontal barrier, 
. Bare soil with vertical barrier, 
. Grass at surface, 
. Grass with horizontal barrier, 
. Grass with vertical barrier, 
. Tree at surface, 
. Tree with vertical barrier, 
. Flower bed at surface, 
10. Flower bed with horizontal barrier, 
11. Flower bed with vertical barrier, 
12. Moisture effect some distance fro 
the barrier/slab edge, 
13. All quantities calculated. 


OMAN ODA AWD = 
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49 
go to 
pg 113 


eo 50 
go to 


45,6 Dg 113 


7,8,12 = 


go to 


9,10 pg 113 
ri t2 


52 
go to 
og 113 





Retrive Values fora 
Particular Case of Post 
Construction Soil Volume 
Change Calculations 


from 
Dg 112 


Initialize And Retrieve Values 


6, = 40 (cm) 


Az 


slab 
i = 0 (soil increment counter) 


z,= 0 (depth to i" increment) 


increment = 5 (cm) 
retrieved, from pg 86 
retrive d,,.., from pg 86 
retrieve HeoilSptvcsary) from pg 110 
retrieve Reap (Zterary) from pg 111 
retrieve h,.4(Zst onsvary) from pg 110 


A(Z st sary a sou 7 San) 


vbar 





53 
go to 
Dg 114 


Initialize And Retrieve Values 


6, = 40 (cm) 
AZ.age = 9 (Cm) 
AZap = 0 (cm) 
i = 0 (soil increment counter) 
z, = 0 (depth to i" increment) 


increment = 5 (cm) 


retrieve d,,., from pg 86 
retrive d,.., from pg 86 
retrieve d..uce from pg 86 
retrieve Ny ee(Z sty cary) from pg 110 
retrieve h,4(Zpteonsyary) from pg 110 


retrieve h,..(Zstyeary) from pg 111 
(Zt yciary = Diree(Z >t verary ) 
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Initialize And Retrieve Values 


6, = 40 (cm) 
AZ, age = 0 (Cm) 
AZ.iap = 9 (cm) 
i = 0 (soil increment counter) 
z,= 0 (depth to i" increment) 


increment = 5 (cm) 
retrieve d,,., from pg 86 
retrive d,,.,, from pg 86 
retrieve h,,..5(Zstyciary) from pg 110 
retrieve h,,..(Ztonsvary) from pg 110 
retrieve h...4(Zpt year) from pg 111 


AZ st ycldry)= grass (Zire vcrdry ) 


Initialize And Retrieve Values 


6, = 40 (cm) 
AZ. age =0 (cm) 
A2Z.ap = 9 (cm) 

1 = 0 (soil increment counter) 
z,= 0 (depth to i'® increment) 
increment = 5 (cm) 
retrieve d,,., from pg 86 
retrive d,.4,from pg 86 
retrieve d,,.,,.e from pg 86 
retrieve hy,(Zst cary) from pg 110 
retrieve h,an(Z»tconsvary) from pg 110 
retrieve h,..(Zstyeary) from pg 111 


D(Z pt ycsary =P (Zptyerary) 





Post construction Soil 
Volume Change 
Calculations With or 
Without Barriers 
(Part | of Il) 


54 ; Yes 
goto —No——<__ is Zz, <= Zm? > 
g115 “ ee 


AH 
AZ cage = AZ cage + INCrement - 7 









h 


a 
V 


) 


114 


es : 
from 
Pg 113 









Convert all Suction in pF 
to centimeters. 
h(cm) = 10°) 


<a 


isi+t 
Zz, = Z,_, + increment 





Y 


1S Ahern)? Neanltibean)? | 





Yes ee oO 
| 
vy 
k, = .33 ace: 
f=.5 f= .8 
e..=center lift from pg 88 e..= edge lift from pg 88 
edge me , d.bar eigen eine = dhbar 











Calculate Horizontal Velocity 








a 








(7), = -neteed 252) 










.)- h - 
V,(z,) = Ko(Zi) Ihe of: ae ‘i 
i? *ve/dry 


Cmcle T duce ) 
Calculate Suction at Foundation edge 
=V, (2,) race 


)-e K,(z;){ho| 

















ee (Z;, lean) = hy (Z;, dry 


i vee he(Z)) = Negge(Zir trary 
<Cis ze 7 
os: ee ee No_ 
he(Z)) = DyaplZitrvary 


AV FLY, 
&2) 7 - (2 Noau( 222) hy(Z;) = Ayan (Zisbrevary 











— 











. ee, 
No _isz,<= 40cm? -Yes 


eee 








Post construction Soil 
Volume Change 
Calculations With or 
Without Barriers 
(Part Il of Il) 


Az AZ at AZ 


total — 


slab edge 





No 


AH 
AZi4, = AZy + Increment - 7 


44) 
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Convert all Suction in pF 
to centimeters. 
h(cm) = 10°F) 






Z,; = Z,, + increment 











h(Z) = hy (Zt He(Z) = ep (Zi teary 


( 


ora 










AV 
“Y) = —7,,(Z,)log,, 


h 






h,(Z;) 
h, (z;) 











APPENDIX B 


SOFTWARE PROGRAM DESIGN LAYOUT 
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Program Design Layout 








Program Design Layout 






































A three dimensional rectangular (X,Y,Z) coordinate system 1s used to define 
foundation layout. The Y direction is defined by some input angle wrt true 
North. The Z direction is depth of soil (positive is down). Every part of the 
foundation structure must be in the positive X,Y plane. The origin is defined 
by Latitude, Longitude, and Elevation. 
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Program Design Layout 





Se ee nn ee ae 
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Program Design Layout 


oy th GRAPHIC | 
: IDISPLA Fi Z 5 ~ 
Foundation Layout | Moisture Effect pu 














AEP Sa 


vyioisture Hit 


Program Design Layout 


fas Gn a5 with horizontal barrier - 
6. Gn 5 with vertical barrier 
?. Tree near foundation edge 


i 9. Tree with vertical barrier 


9. Flower bed near foundation edge 


10. Flower bed vath horizontal barrier y il 


nT. ‘Flower b bed with vertical barrier _ 


12. Measured suction profiles 


orizontal 
Barrier 
Tee wi roots 
under slab. 
Leaking 
Pipe Vertical 
Barrie 


Live Oak _——+-* 
Tree Removed 
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Program Design Layout 





Number Initial Final : Surface 
of Suction Suction Suction - Conditions 
Measurements' Profile(pF)) Profile(p F) (Vegetation/Bare) 


Case 
Number 


AP Ah et ANAL ARMA At LORNA NAAN RS te bd 


Nn a 


Vegetation 
Vegetation 


OR TORR TR eR ~ ; 


| Vegetation 


ekeths iaribethithorsthett ith elteatheteetieteees 
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Program Design Layout 





Suction _ Depth to Suction 
Measurement ‘Measurement Measurement 


Number 7 (ft) | (pF) 


3.0 
12.0 


we rene Are 


| 


fon G0 ND 


"16.0 


i 


2s 
PO a lO 
Final 
Final 
Final 


ei ND 


Initial 

_ Initial 

Initial 
Initial 
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Program Design Layout 
GRAPHIC D 


0 a SS SSS ee 
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Program Design Layout 





oe a 


| GRAPHIC DISPLAY 


4 ar £ 
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Program Design Layout 
CCS ae ; eran 


= ieee 


[ola] 


Origin 


Goto Soil Layer 
Property Input —_ = 


eee on Soil ae cs al erties 


%LL “PL %#200' %2zp “PI ‘Fine Slaw CEC | 


0.7104 40.6806 0.163) 0. 190650 120, 
0.6789 _ 0.7682 a 0.062426 120 


0.000681 2.96E-06 
0.000648 2.81E-06. 
0.000621 2. 70E-06 


_| _0.000635._2. T6ED6 








| 15Grasswihorizbar — 


1o0G wi vert, | A 
a Wh Vert, 
AAAS AC ASIAN EEE AN ERE EEE ERR SIARNNS 


LH es ay, 


Grass with 4ft Vertical Barrier 
Section A-A 


Suction Profile At Slab Edge With Barrier. 
Suction Profile Under Slab in June. 
Suction Profile At Time of Construction. 
———__—+— Equilibrium Suction Profile 3.31 pF. 
— — Maximum Suction Profile. 
SS Minimum Suction Profile. 
Gvbar = 4ft Depth of Barrier 
_m = 10.3ft Depth to Constant Suction. 
Mm = 10.5ft Edge Moisture Variation Distance 


aoe Pl ett ot ae 


: 
| 
| 
| 
| 
| 





1, 
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Poway Ss Ta LIS Le 


1-5 Grass wi horiz.bar. 


vrahaa: Sn En RE em mmo 


August _j])}orsets 
Slab Constructed in | January 
Driest Month is [__sune__| 


b Grass wi vert. bar 


_1-12 Measured 


Section A-A 


Slab at time 


of construction Soil at time 


of construction 


Soil movement 
through the month 
of August. 
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Program Design Layout 





Royo eb aa ky 


a a ae 


Section A-A $12 Mea sured 13 


‘ 
ss 


Suction (pF) 


weedy 


j 
| 
l 
5 
5 @ 
l 3 
| 210 = 
i =~, 
O 
i a5 
| c= 
i 
i 
J 
\ 
| 
1 


———= - Suction profile due to the flower bed. 
----—- - Equilibrium suction profile. 


anne 8 nna ne 


Diener Mea RI ANAND ALAMOS SAAD ih ath Sh VIDNDAG HOON ATOR ARDY MAG NR SAA SANA A MAIEV ION NL A I SRSA IODA NO AAMAS Ee 


Se S Fo 3 : 5 Seid An ecco Ores ate. ES 


Cm = 6.5 (ft) Edge moisture variation distance. 
Zim = 13.8 (ft) Depth to constant suction. 
hin = 5.57 (ft) Equilibrium suction at the surface. 


d = 4.0 (ft) Depth of the flowerbed. 


zone 
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Program Design Layout 





| GRAPHIC DISPLAY = 





APPENDIX C 


SITE PLAN OF MOISTURE EFFECT CASES 
USED IN SAMPLE CALCULATIONS 
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OAV DOT PS Yi stebt hag pe 


LLP AES AY PIAS AOL iO ie 


SPRUE RRC EUR EERIE, ISO BON ARNT ETE D OM NEALE ELIE GNG SRST SLE RO OO E LI EI OR NE LEI ND RINE TN IL ETE NL IONE, EIEIO TIA ROP LIN LILES LILLIE GGG 
Be 


GRAPHIC DISPLAY 





Depth of Tree 
Root Fiber 


is 20° 
on i 
35' 
Qak 
Tree 


20° 


Depth of 
Flower Bed 
Zone is 4 


rhe tai ntebre tentetnlenncnacn. na pettiness he ihe nnn heeenabainnnnrpinfanme ttn nnn rene nite mbt anhA nN At ntshtehttnidant te Bt BD NN Nat DO ANN AO OOD tee Nein Hanon neti 
Nabe ones ve sits. =o ~ . <= a 


Program Design Layout 





Foundation Example 
Used to Illustrate Program 
Design Layout 4 0 


65" 


A* Sidewalk 


Point of 
i Water Source 


Vertics)] Barrer 
(Length of East Edge) 








{20° 


Live ae 


Tree Remoaved 


Grassy Lawn 
{Around the Entire Building} 


‘ 
t 
i 
+ 


-> 


' ee 
cree an 


1 


Water Source 





Flower Bed 


Project Name: Sample 

Project Date: 4 April 1997 

Project Engineer: Donald D. Naiser, Jr. 

Site Location Address: 400 Marion Pugh #1905 
College Station, TX 77840 

Approximate Location: Lat. 29 15° N 
Long. 96 28° W 


eA Dey ty ed Napanee pareraapaeren one Reoperqwovenscoay: i eae x ~ gn mera Se Sl Ie as ae a ee eee ae 


tne ant none nbint tebe 


GSO 


8 LEE AI CIP ALR AI, AREAL AE 
oer 


32 





APPENDIX D 


REQUEST FOR PERMISSION TO USE COPYRIGHTED MATERIAL 
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Donald D. Naiser, Jr. 
400 Marion Pugh #1905 
College Station, Texas 77840 


August 27, 1997 


Permissions Department 
John Wiley & Sons, Inc. 
605 Third Avenue 

New York, NY 10158-0012 


REQUEST TO USE FIGURES COPIED FROM SOIL MECHANICS 1969 


Iam a master of science candidate at Texas A&M University and I will be finalizing my 
thesis in the next two weeks. Our library sends our theses to University Microfilms Inc., for 
preparation of a microfilm copy of the document. As I am sure you are aware, UMI retains the 
right to sell both hard and microfilm copies of the document. 


My reason for writing is to request permission to use the figures listed below exactly as 
they appear in SOJL MECHANICS, by Lambe and Whitman, copyrighted by John Wiley & 
Sons, Inc., 1969. With your permission, I will use the following figures in my thesis: 


Figure 4.7 “The structure of serpentine.”’ Found on page 47. 
Figure 4.8 “The structure of kaolinite.” Found on page 47. 
Figure 4.9 “The structure of pyrophyllite.” Found on page 47. 
Figure 4.10 “The structure of muscovite.” Found on page 47. 
Figure 4.11 “Sheet silicate minerals.” Found on page 48-49. 
Figure 5.9 “Soil particles with water and ions.” Found on page 55. 


oo! ae 


One chapter in my thesis describes the magnitude of the problem of designing 
foundations built on expansive soils. The figures described above will be used to illustrate the 
molecular structure of fine grained soil particles. The thesis has been accepted and the author 
proofs indicate a 1997 publication date. The thesis title is PROCEDURES TO PREDICT 
VERTICAL DIFFERENTIAL SOIL MOVEMENT FOR EXPANSIVE SOILS. 


Please inform me of any necessary action to prevent future publication complications. | 
can be contacted by phone at home (409) 543-8412 or at work (409) 845-9919. Thank you very 
much for your effort in resolving these issues. 

Sincerely, 


Ope LOGE 


Donald D. Naiser, Jr. 
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VITA 


Donald David Naiser, Jr. is a citizen and resident of the United States of 
America, whose permanent mailing address 1s, 
Rt. 3 Box 203-C 
El Campo, Texas 77437. 
He obtained his Bachelor of Science degree in Civil Engineering, from Texas 
A&M University in 1991. 
Since graduating from Texas A&M University, he has been commissioned as a 
United States Naval Officer as a member of the Civil Engineer Corps. His time served in 
the Navy has been exclusively with the SEABEES performing contingency construction 
training as well as constructing naval facilities. His tour highlights have been as the 
Officer In Charge of the Elevated Causeway System (ELCAS) in Amphibious 
Construction Battalion ONE and as the Officer In Charge of Construction Battalion Unit 
418. 
In 1996 he was accepted at Texas A&M University to pursue studies toward the 
degree of M.S. within the Geotechnical branch of Civil Engineering. 
In 1997 he became registered as a professional engineer in the state of 


Washington. 


“Blessed are the people whose leaders can look destiny in the eye without 
flinching but also without attempting to play God.”, Henry Kissinger. His career 


interests include achieving excellence through leadership in any Christian form. 
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